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CHAPTER  I 
INTRODUCTION 

A  production  system  raay  be.  dofinea  os    an  orderly  collectioa  of 
elements  of  production  v;hich  are  ccmbined  and  operated  in  a  certain 
desired  fasbioa  co  produce  goods  of  economic  use  or  to  result  into 
an  efficient  service.   Tlie  basic  elements  of  a  production  systeir.  are 
man,  machine  and  materials. 

The  production  system  can  be  analyzed  in  terms  of  the  nature  of 
flow  and  the  manner  th.e  jobs  are  handled  in  shop.   The  two  broad  categories 
into  which  a  production  system  can  be  classified  are  batch  product  ion 
and  contiLvacus  production  systems.   In  a  batch  production  system,  ithe 
jobs  each  of  which  consists  of  a  certain  lot  size  are  processed  oa 
the  various  machines.   A  job  shop  manufacturing  electric  motors  is  an 
cx9japle  of  batc]\  production  system.   The  distinguishing  feature  of 
such  systeri  is  tb.e  multi-purpose  nature  of  the  machixies  vised.   All 
machines  in  the  sh.op  are  universal  and  no  machine  is  used  for  the 
specific  purpose  of  a  particular  job.   These  machines  nay  be  arranged 
according  to  their  types  and  each  job  may  have  different  routing. 
Therefore,  such  a  system  is  referred  to  as  shop  production  system. 
In  continuous  production  system,  highly  specialized  nature  of  the 
machines  chatacterizes  such  system.   Car  maauf acturing  industry  is 
an  example  of  continuous  production  system.  T\\e   machines  perform  a 
specific  operation  on  each  product  which  form  a  line  and  approach 
individually  towards  a  machine.  Also    the  Tnachines  are  arranged  in  a 


parLicular  order  to  coaforra  to  the  order  of  operations  to  be  performed 
on  Lhe  jcl).   Such  a   aiachinc  arrangement  ±s    commoaly  called  an  assembly- 
line  and  the  system  is  referred  to  as  line  production  system. 
■^'  -^   Job-Shop  SchedullBg 

Tliis  thesis  is  concerned  with  the  shop  production  system.   In 
such  a  system,  the  shop  consists  of  a  set  or  group  of  working  centers 
and  a  set  of  operations  to  be  performed  in  a  specified  technological 
ordering  at  one  or  more  working  centers.   Each  center  has  a  finite  number 
of  machines  which  are  identical,  as  far  as  the  nature  of  work  demanded 
at  each  working  center  is  concerned.   The  complex  behavior  of  a  job 
sKop  can  be  described  as  follows: 

After  a  customer  order  is  received,  the  engineering  department 
decides  the  nature  of  the  work  to  be  performed  on  different  machii'ies 
and  specifies  the  technological  ordering  of  machines  to  perform  the 
work  on  the  o.Jer,  hereafter  referred  to  as  a  job.   Engineering 
departmieat  also  sets  the  time  limit  by  which  each  job  has  to  be 
completed,  kno\-m  as  the  due-date  of  a  job.   In  some  cases,  a  job  may 
have  other  attributes  of  interest  such  as  monetcry  value  or  it  may 
have  a  priority  because  of  its  regular  potential  customer.   The 
information  necessary  for  performing  operations  of  a  job  is  sumiuari^ed 
on  a  paper  generally  tenned  as  a  job  file  and  accompanies  each  job 
throught  the  shop.   After  a  job  has  arrived  in  the  shop  the  ma:;-.ine 
on  which  the  first  operation  is  to  be  processed  according  to  the 
corresponding  routing,  is  checked  for  idleness.   If  that  riachine  is 
processing  another  job,  then  the  newly  arrived  job  joins  a  q,(:euo  before 
the   work  center.   However,  if  trie  macliine  is  idle  and  rhc  infor'aation  of 


the  next  job  to  arrive  is  absent,  the  madiine  starts  to  prf^pare  for 
the.  operation  while  the  job  waits  at  the  work  center.   The  tinie  spent 
in  the  preparation  of  the  aiachine  for  processing  that  job  is  knovn  as 
set~up  tifne  of  a  job.   After  set-up  tir.ie  is  over,  a  r.iachine  processes  the 
j;ob  for  a  pre-specified  operation.   The  time  interval  between  end  of  the 
set-up  and  the  corapletion  of  processing  the  job  is  called  operation 
ptrocessing  time.   The  job  is  then  transported  to  another  machine 
specified  in  the  job  file  for  performing  its  second  operation.   After 
all  operations  of  a  job  are  completed,  it  departs  from  the  shop  as  a 
CDHipleted  job. 

Siich  a  complex  behavior  of  a  job  shop  necessitates  assumptions 
of  some  theoretical  restrictions  if  experimental  work  is  to  be  carried 
out  on  job-shop  models.   A  simplified  job-shop  model  is  considered  to 
have  the  follov/ing  assumptions  [2,  12,  23,  29,  ^'C]: 
r,.  Assumptions  related  to  jobs 

1,1'  The  arrival  pattern  of  jobs  in  the  shop  is  probabilistic  in 
nature.   The  governing  probability  distribution  remains  un- 
changed throught  the  length  of  an  experiment  and  is  not  c^ffected 
by  shop  conditions. 
1.2^  The  arrival  of  a  job  to  the  shop  includes  the  infom-.aticn 

necessary  for  processing  such  as  processing  titnes,  job  routing 
and  due-date  are  kno-.vn. 
r. 3  Even  though  a  job  may  represent  a  lot  consisting  of  several 

individ^ial  parts,  no  lot  is  processed  op.  more  than  one  nachine 
at  a  time,  i.e.,  n-j  lot  splitting  i.s  permitted,  or .. .overlapping 
of  operations. 


1.4  Job  routing  is  fixed  and  no  alrernate  routings  nre  permissible. 

1.5  Each  job,  once  started,  is  procesi^ed  untiJ  its  coriipj.etiou; 
i.e.,  no  order  cancellation  occurs. 

1.6  Each  job  may  wait  between  machines;  i.e.,  in-process  inventory 
is  allowed, 

1.7  No  interaction  is  allowed  among  jobs  at  any  work  center.   Jobs 
of  nearly  identical  set-ups  do  i:ot  accumulate  in  order  to 
save  set-up  times  at  given  work  center. 

2.  Assumptions  related  to  productive  elements  (man  and  machine) 

2, J   Each  machine  center  has  a  single  queue  of  jobs  to  be  processed. 

2.2  No  machine  may  process  more  than  one  job  at  a  time. 

2.3  Each  machine  is  continuously  operating  in  time,  i.e.,  time  is 
not  divided  into  shifts,  days  or  weeks. 

2.4  Machine  breakdowns,  manufacturing  errors,  employee  absenteeism, 
material  shortage,  design  changes  and  such  interruptions  are 
not  permitted. 

2.5  The  system  is  either  labor  limited  or  machine  limited. 

3,  Assumptions  related  to  operations 

3.1  An  operation,  once  started,  is  perfotxaed  to  its  completion;  i,e., 
no  pre-emptive  priorities  are  allowed. 

3.2  Each  operation  can  be  performed  by  only  one  machine  at  a  tirae. 

3.3  Operation  processing  times  are  independent  of  the  order  in 
which  the  operation  are  performed.   In  particular,  set-up  times, 
if  any  J  are,  are  included  in  processing  times  and  are  sequence 
independent.., 


3.4  Each  operation  of  a  job  inust  be  coinpletcd  before  Its  succeeding 
cperation  can  begin  i.e.,  no  phase -l.-ipping  is  allowed. 

3.5  Transportation  times  from  one  machine  to  another  are  not 
allov/cd  or  inclaJed  in  processing  times. 

Inspite  of  the  above  simplifications  it  is  presumed  that  the  model  still 
rata.ins  the  essence  of  a  realistic  job-shop  which  is  much  moie  complex 
in  nature  than  the  one  v/ith  assumptions. 

Basically,  the  job-shop  model  may  be  viewed  as  a  queueing  system 
involving  a  network  of  t\'70  or  more  distinct  work  centers  (servers). 
Each  job  is  required  to  be  processed  by  a  finite  number  of  machines 
according  to  the  specified  routing.   The  objective  is  to  schedule 
the  jobs  such  that  a  desired  criterion  is  optimized. 

Tlie  tenn  industrial  scheduling  is  used  to  refer  to  the  planning 
of  a  Ij.mited  manufacturing  capacity  to  achieve  certain  goal  over  a 
par:Iod  of  time.   Such  a  procedure  of  planning,  a  schedule  is  necessary 
in  order  to  utilize  the  limited  production  capacity  effectively. 
Scheduling  can  be  classified  into  two  categories  depending  on  the  time 
horizon  used;  namely  loiig-tcrm  and  short-term  scheduling.   Flanning 
the  gross  number  of  worlanen  and  number  and  type  of  machines  required 
to  manufacture,  say,  two  turbo  generators  over  a  period  of  one  year 
is  an  example  of  long-term  scheduling.   While  planning  a  number  of 
operators  in  a  job  shop  on  a  day-to-day  basis  is  an  example  of  short- 
term  scheduling.   This  type  of  scheduling  form.s  the  main  topic  of 
discussion  for  this  thesis. 

It  is  noted  that  the  terms  scheduling,  dispatching  and  sequencing 
have  bee:;  used  i  nterclringeably  by  many  authors.   Each  of  the  above  terms 


is  defined  as  considered  in  this  Lliesis.   Scheduling  is  f.h.e  l.iiuing  of 
the  operations  performed  on  each  machine  so  as  to  optimize  some  desired 
criteria.   However,  dispatching  is  the  seqj.ieuciug  of  jobs  on  a  facility 
without  taking  into  consideration  the  effect  of  tln^.e.   Dispatching  is 
also  called  sequencing. 

In  order  to  run  a  shop  effectively,  some  scheduli.ng  rules  are 
required  instead  of  the  rules  of  thumb  which  have  been  used  for  years 
in  actual  shops.   Scheduling  rules  help  in  selecting  a  job  from  a 
number  of  jobs  v/aiting  to  be  processed  on  a  machine.   These  rules  are 
also  referred  to  as  decision,  dispatching,  sequencing,  loading  or 
priority  rules.   Some  of  these  rules  are  called  heuristics.   Scheduling 
rules  are  laeans  of  dispatching  procedures  to  generate  a  schedule. 

Scheduling  of  jobs  in  a  single -channel  queueing  system  has  been 
investigated  thoroughly  by  analytical  methods  but  unfortunately  no 
analytical  solutions  are  reported  for  network  configuration  systems. 

A  number  of  investigators  have  therefore  resorted  to  experimental 
investigations  by  digital  computer  simulation.   Although  simulation 
does  not  provide  vzith  the  optiviial  solution,  it  provides  some  information 
v/hich  could  motivate  analytical  work  on  network  systems.   Considerable 
amount  of  research  along  this  line  has  been  reported  in  a  summary  fashion 
in  Chapter  II. 
1 . 2  Proposed  Research 

The  purpose  of  this  research  is  four-fold. 

1.   To  provide  a  complete  aitd  comprehensive  reviev/  of  the  most 
of  the  research  on  job-shop  simulation.   An  effort  has  been 
made  to  foi'm  a  tabulated  summary  of  simulation  expcrimeni.s 
including  shop  characteristics,  experimental  conditions, 
scheduling  rules,  measures  of  performance  and  results  obtained. 


2.  To  develop  a  coinputational  procedure  for  job  Uue--date  and 
evaluate  it:  in  coiupacison  to  the  best  procedure  proposed 
by  Conway  [12]  . 

3.  To  test  the  computational  ef;ficLeacy  of  GASP-'II  in  handling 
simulation  experiments  for  job-shop  situatlous. 

4.  To  study  the  effect  of  the  change  in  the  processing  time  dis- 
tribution from  exponel^tial  to  Erlang  with  different 
parameters.   The  rules  tested  are  those  v/hicli  have  been  found 
by  Conway  [12],  in  general,  superior  than  others,  namely, 
shortest  imminent  operation  time  ['RUL10  5]  ,  remaining  job 
slack  time  [RUL114]  ,  and  reinaining  job  slack  to  remaining 
number  of  operation  [RUL226]. 


CHAPTER  II 

LITEMTURE  REVIEW 

This  chapter  is  devoted  to  a  comperhensivti  revievz  and  siuTtmary  of 
most  of  the  research  done  so  far  on  the  performance  of  the  Vcirious 
scheduJ-lng  rules  applicable  to  job  -shop  sch.eduling. 

^- '  1   Scheduling  Rules 

Scheduling  rule  may  be  defined  as  the  procedure  by  v/hich  a  job  is 
selected  from  a  queue  for  scb.Gduling.   Although  the  term  "scheduling 
rule"  is  used  in  this  paper,  the  rules  may  be  referred  to  as  sequencing, 
loading,  dispatching,  priority,  heuristic,  on  decision  rules.   Scheduling 
rules  can  be  classified  into  various  classes  depending  on  the  amount  and 
type  of  inf ormiition  required.   Four  classes  into  i.-'hich  the  rules  can  be 
grouped,  are  as  follows  [46] : 

1,   Local  Rule;   A  local  scheduliiig  rule  utilizes  only  the  information 
associated  v^ith  the  jobs  in  the  queue  from  which  a  job  is  to 
be  scheduled  next.   'llie  status  of  other  m.achines  i.s  not  con- 
sidered while  the  selection  is  made, 
^'   Global  Rul'.^-;   While  selecting  a  job  from  the  queue,  a  glc'bal 
scheduling  rule  considers  the  information  avciilable  in  the 
queue  as  well  as  the  status  of  other  machines  and  also  the 
attribtites  of  jobs  waiting  in  the  various  queues  in  the  shop. 
■^  •   Static  Rule;   When  the  assigned  priority  index  does  not  change 
with  time,  the  rule  is  said  to  be  static.   For  ~.xa:nple,  first- 
come  first-served  and  shortest  ope'^'atjon  prccessirtg  ti;r>c  rules 
are  scatlc. 


^'      Dyriamic  Rule;   When  the  priority  of  jobs  waiting  in  the 

qurue  changes  with  tine,  the  rule  is  referred  to  as  a  dynamic 
rule.   For  example,  the  scheduljng  rule  in   which  prjority 
is  given  inversely  proportional  to  tlie  nuraber  of  jobs  wiiiting 
in  the  queue  of  the  subsequent  machine,  according  to  machine 
ordering,  is  a  dynamic  rule. 

Various  scheduling  rules  (70  rules)  which  came  to  our 
attention  xvhile  surveying  the  previous  work,  are   listed  in 
Appendix  A.   These  rules  have  been  carefully  reviewed  and 
defined  including  mathem^itical  expressions  except  heuristics. 

2 . 2  Measures  of  Performance. 

Measures  of  performance  are  the  basic  criteria  upon  v;hich  the  test 
of  goodness  of  the  rules  is  based.  Two   important  categories  into  '■-•bich 
thiese  can  be  classified  are  [12]:   (1)  those  which  are  related  to  rhe 
jobs;  and  (2)  those  which  are  related  to  the  shop.   However,  \/hen  a 
measure  of  p-irformance  is  a  function  of  cost,  it  m.ay  be  related  to 
attributes  ot  both  shop  and  jobs  together.   Thie  summ.ary  of  the  various 
m.easures  of  performance  which  have  been  most  frequently  considered  in 
job-shop  scheduling  research  is  given  belov/.   For  convenience,  the 
notation  used  are  summarized  below 

n     sample  size 

D,  due-date   of    iob  j 

C.  ccmol-Ction    time    of    iob    i 

J  - 

L.  lateness   of   job   j 

L  ineau  lateness   of    late  jobs 


JO 


V^  variance  of   lateness 

Li 

N      number  of  late  jobs 

T.      tardiness  ot  job  i 
J 

T      conditiona].  tardiness 
c 

T  mean  tardiness  of  jobs 

A.  delay  factor  of  job  j 

W  mean  job  waiting  time 

F.  shop  flow  tirae  of  job  j 

F  mean  of  shop  flow  times 

V  variance  of  shop  flow  times 

M.  (t)  aum.bcr  of  machines  busy  at  time  t 

N(l)  total  number  of  jobs  waiting  in  queues  at  tirr.e  t 

N(t)  total  number  of  jobs  in  shop  at  time  t 

Q(c)  total  shop  work  at  time  t 

Q'''(t)  total  iimuinent  operation  work  at  tim.e  t 

Q(t)  total  shop  v7ork  uot  yet  started  at  time  t 

I .   J cb  Charac teris  t i cs 

Job  cbaracterisrics  are  measures  of  performance  which  give  informatioa 
about  cattributes  of  jobs  such  as  latraess,  waiting  tune  in  'lueueSj,  ai>d 
shop  fJ.ow  tii~e.   The  following  16  measures  of  perfoniiance,  found 
mostly  in  [12,  23,  23],  are  given  be  lev/. 

MESlOl;   JOB  LATENESS 

Job  lateness  is  defined  as  algebraic  difference  betweexa  job 

completion  tiioe  and  job  due-date. 


C.  -  D, 
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MES102;   JOB  X'MDINESS 

Job  tardiiicss  is  defined  as  the  positive  lateness 


T.  -  mox  \0,    C.  -  D.] 
J  J    3 


-  max    [Oj    L. ] 
IVLES103:      MEMT  LATENESS   OF  JOBS 


L  =      I     L./n 


MES104:   VARIM'CE  OF  LATENESS  OF  JOBS 

n   (L  )2 

j=l 

MES105:   CONDITIONAL  TAPJDINESS 

Conditional  tardiness  is  defined  as  mean  tardiness  of  all  tardy 

(laue)  jobs, 

N 
T   =   y   T./N 

MES106:   MEAN  TARDINESS  OF  JOBS 

N         NT 


L 


y   T./n---^ 


J-1 

MES1G7:      TOTAL   iNTJRBER  OF  JOBS   LATE 
■■!otal  nivuber  of   jobs   for  \7hich    lateness   is   positive. 


N  =      y     X. 


vjhere 

X.    =   1,        if        L.    >      0 
J  '  1 


-  0,        if        L,    <      0 
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MES108:   SHOP  FLOW  TIME  OF  JOB 

Shop  flew  time  of  a  job  Is  defined  as  the  difference  between  the 

job  compleliloa  time  and  the  job  release  tlae. 


F.  -  C.  -  r., 
J     J     J  ^ 


MES109;   M£M  SHOP  FLOW  TIME  OF  JOBS 

n 

F  -      I      F  /n 

>IES110:      VAEIANCE   OF  SHOP   FLOW  TIME   OF  JOBS 

V^   =      )      F^/n   -    (F)^ 

>;IES1.U:   DELAY  FACTOR 

Delay  faction  is  defined  as  the  ratio  of  the  shop  flow  time  of  a 

job  to  the  sum  of  its  operations  processing  times. 

K. 
A.  =  (C.  -  r.  J   )    t., 

MES112:   PENALTY  COST  FUNCTION 

Penalty  cost  function  is  defined  as  weighted  sum  of  the  mean  job 

lateness  and  the  standard  deviation  of  its  lateness. 


Z  -  aL  -I-  b/v" 
where  a,  b  cost  parameters  (constant). 
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MES113;   PENALTY  COST  FUNCTION  BASED  JOB  LATENESS 

Pan.-iJ.ty  cost  fuxiction  based  on  job  tardine<is  is  defliied  as  the 

weighted  job  lateness. 


P  =  c  (C.  -  D.) 
J     J 

where  c      weighting  factor  such  that 


'c^,,   if  D.  >  C. 
t-        2  3 


c,  ,      if   D.  <  C. 
L        J    J 


MES114;   MEAN  JOB  WAITING  TIME 

Mean  job  v/aiting  time  is  defined  as  the  time  elapsed  from  the  raomerit 

job  arrives  until  its  processing  on  a  machine  starts,   Tliis  neasuce 

of  performance  is  based  on  all  jobs  completed  during  the  simulation 

run. 

K. 

j=l  k-l        ^'- 
MES115:   COST  FUNCTION  BASED  ON  TARDINESS 

Cost  function  based  on  tardiness  is  defined  as  the  weighted  sum  of 
tardiness  and  the  square  of  tardiness  of  all  jobs. 

n  „ 

where 

X^ ,  K^  and  K„  are  constants. 


MES207;      TOTAL  OUTSTAI^DING    (REit/vINING)    SHOP  WORK. 

Total   reraaiiiiug   shop  work   is   defiued   as    Llie  sum  of   processing    tiirios 

oi"    all   operaLions  which  have  not  yet   started   for   all  jobs   in   the 

shop   at    time   t. 

K. 

jesh.op  k~k 
MES208:   MEAN  TOTAL  OUTSTANDING  SHOP  WOEK,  Q 
MES209;   VARIANCE  OE  TOTAL  OUTSTANDING  SHOP  WORK  V- 
MES210;   IMMINENT  OPERATION  WORK 

Imminent  operation  work  is  defined  as  the  sum  of  imiainent  operation 
processing  times  for  all  jobs  in  queue  at  time  t. 

QMt)  -'^   I  t  i; 

jGshop  -^ 

MES21.1;   MEA.N  OF  IMMINENT  OPERATION  WORK,  Q* 

MES212:   VARIANCE  OF  IMI4INENT  OPERATION  WORK  V  ^. 

MES213:   TOTi\L  NUl-IBER  OF  JOBS  IN  SHOP 

Total  number  of  jobs  in  shop  is  defined  us    the  sum  of  all  jobs 
v;aitiug  in  queues  and  jobs  being  processed  on  machines  (busy 
machine) . 

N(t)  -  N(t)  +  M,  (t) 

D 

MES_2_J^:      MEAN  NUMBER  OF  JOBS    IN   SHOP,    N 

MiZ3215:      VARIMCE   OF  NU14BER  OF  JOBS   IN  SHOP,   V;; 

MES216;      x^LA.CHINE   ICLE  TIME 

Machine  idle   tliae   is   defined  as    the  machine   time  betv/een   tlie 

completion  of   an  operation  and   the   arrival  of    the   next   job   for 

processing   cae  of   its    operation.      It   is   m.easured   as   perceatagc 

of.    total  simulation   time. 


-1.:) 


MESil6;  SUM  OF  TAI1531NESS 

When   the  values   of   constants    In   tba   above  neasure  of   perf orr.iance 

(MES1I5)  are  such   tiiat   K     -  K     =   0   and  K^    =--   1,    the  above   cost 

function  becomes    the  sura  of    tardiness   of   ail   late   jobs   such    that, 

n 
Z  =      y      T. 

^^-      Shop  Characteristics 

Shop  characteristics  are  measures  of  performances  vhlch  give  information 
about  conditions  of  the  shop  at  any  glveii  time. 
MES201;   TOTAL  NUIvBER  OF  JOBS  IN  QUEUES 

M 

N(t)  =  y  N.,  (t) 

m-1      J'^ 
MES^_2_:      I'EAN  NIHBER  OF  JOBS,   N 
■^■ES203:      VAEIANCE   OF  NUIIBER  OF   JOBS,    V 
MES204:      TOTAL  SHOP  WORK 

Total   shop  fvork  is   defined  as    the  sum  of  processing    times    of    all 
operations   of  jobs   in    the  shop   at    time   t. 

K. 

Q(t)    =  if      »:;, 

jeshop  fe=l       -^ 


KES205;      MllAN  SIiOP  WORK,   Q 
tIESZOo:      VARl^-XE    OE   SHOP   WORK,    V 


Q 


I'fESZO?:   MCHT.NK  UTILIZATIONS 

xMachines  utilizatioa  is  defined  as  tlie  r/. Lio  of  working  time  of 

machines  to  the  total  available  time.. 

2,3  ^i mulation  Experiments 

Considerable  experimentations  for  testing  the  different  schedulLng 
rules  have  been  conducted  by  several  investigators.   This  section 
summarizes  the  experiments  which  are  worth  while  foe  consideration 
before  any  new  experiment  should  be  designed. 

The  experiments  consists  nainly  of  two  different  natures;  (1)  machine 
limited  systems;  and  (2)  .labor  limited  systems.   Comparatively  very  little 
work  has  been  done  on  labor  limited  systems.   Rowe,  Allen  and,  LeCrande 
[53,  1,  41]  are  those  who,  to  the  authors  knowledge,  have  investigated 
such  systems.   Incidently,  the  experim.ents  have  been  conducted  using 
actual  shop  settings  of  General  Electric  Company,  and  Hughes  Aircraft 
Company.   On  the  other  h.'md,  in  the  field  of  machine  limited  systems, 
Conway  and  his  associates  at  Cornell  University,  have  carried  out  several 
simulation  experiments  to  study  the  effects  of  scheduling  rules  in  Job 
shop.   An  extensive  si'.udy  has  been  conducted  by  Conway  [12]  at  Rand 
Corporation.   In  this  work  9  2  scheduling  rules  has  been  tested  in  a 
maintainance  shop  of  United  States  Air  Force.   Following  is  a  list  of 
16  experiments  which  were  carefully  reviewed.   The  summary  of  each 
experiment  together  v/ith  results  obtained  are  described  hereunder. 
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Experiment  Invei^tlgator  Reference 

Number 

1  Baker  and  Dzielinski  (1960)  [4] 

2  Conway,  Johnson  and  Max^/oll  (1960)  1'17] 

3  Kurtani  and  Nelson  (1960)  [^0] 
h  Conway  and  Max\^^ell  (1982)  iii3] 

5  Gere  (1962)  [29] 

6  LcCraade  (1963)  [41] 

7  Com.7ay  (196A,  63)  [12,16] 
3  Carrol  (1965)  [11] 

9  Or kin  (1966)  [52] 

10  Old '/.ley  (1966)  [lUj 

11  Eilca  and  Hodgson  (1967)  [23] 

12  Hoodie  and  Roberts  (1967)  [45] 

13  Nelmcler  (1967)  [50] 
1.'4  Thompson  (1967)  160] 
15  Eilon  and  Cotteril  (1968)  [22  1 


EXVERIMF.NT  1 
INVESTIGATOR         Baker  and  Dzlelinski  [4] 
COMPUTER/PROG.  LANG  IBM  70A/- 

I   SPxOP  CHARACTERISTICS 

1.  Shop  structure 

2.  Shop  type 

3.  Number  of  groups 

4.  NuTober  of  ma  chines /group 

5.  Shop  utilizatioa 
II   EXPERIMENTAL  CONDITIONS 

1.  Dlstr.  01  job  arrivals 

2.  Distr.  of  job  routings 


Network  of  single  queues 

Job  shop 

9-30 

1 
80% 


3,   Distr.  of  Number  of 
operations /job . 


Uniform 

Normal,  c  =  average  no.  of  operations 
per  job  is  2  to  10. 


4.  Distr.  of  processing  time   Negative  exponential,  1/c  -  aiean 

5.  Number  of  runs  — 

6.  Initial  Shop  State  Empty 


7.  Initial  nuniber  of 
released  jobs 

8.  Run-in  length 


10  and  20  unit  loads  (A  unit  load: 
aggregate  of  jobs  involved  requires 
one  unit  of  expected  processing  time 
from  each  facility 


9.   Number  of  jobs  considered 
for  collecting  statistics 

10.  Run-out   length 

11.  Jobs /run  or   simulatioa 
tivie 


5  unit    loads 


12.      Due-date  Biultiplier 
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EXPERIMENT  1  (conr.) 
Ill   Schaduliug  rules:   RULlOl,  102,  105,  109,  111,  120,  122,  123 
IV  N.easuroa  of  p;^rfoi'naac:e:   MKS109 
V  Results: 

1.  Job's  total  expected  processing  tirne  and  the  nuiiiber  of 
procesalng  operations  per  job  had  a  significant  effect  at 
5%  confidence  level  upon  the  job's  total  manuf acturi.ng 
tir.-ie. 

2.  If  the  average  of  the  job's  manufacturing  time  is  used  as 
a  iDeasure  of  effectiveness,  tlie  simple  shortest  irmninent 
processing  tirae  rule,  is  the  best. 
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INVESTIGATOR 


EXPERIMENT  2 
Conway,  Johnson  and  Maxwell  [17] 


COMPUTFR/PROG.  I:\NG.  IBM  650/- 

I   SHOP  CiiARACr;uliSTICS 

1.  Shop   ;3trucl:ure 

2.  Sliop   f.ype 

3.  Number  of  groups 
h.      Number  of  machines /group 
5.   Shop  utilization 


Network  of  single  queues 
Job  shop 

5 

1 

Light,  liiediuin  and  heavy  load  v.ith 
average  total  processing  time  of 
60,  120,  and  200  time  units 
respectively. 


II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals 

2.  Distr.  of  job  routings 


3.   Distr.  of  Number  of 
operations /job. 


Expected  number  of  operations  on 
each  machine  is  equal. 

Geometric,  p  -  ->-,  7  operations 
maximum. 


A.   Di;jtr.  of  processing  time   Total  allowable  processing  time 

per  job  is  normal 


5.  Number  of  runs 

6.  Initial  Shop  State 

7.  Initial  number  of 
released  jobs 

8 .  Ru  n-  in   1  eiie  t  h 


Non-empty,    depending   on  shop   load 
FuiTction   of   shop    load 


9.      Number   of   jobs    coiisidered      1.      First   35    completed  jobs 

for   collecting   statistics      2.      First    85    completed  jobs  -1-   last 

15  jobs   forced   out 
3  .      A  J.  .1    1 J  0  j  ob  3    c  omp  1  e  t  ed  . 


10.       Riin-out    length 


15   jobs 


11,      Jobs/run   or   simulation  100   job; 

time 


12,      Due-date  multi.nlier 


EXPERIMENT  1      (cont.) 

Ill   Scheduling  rules:   RULlOl,  109,  111,  119,  V'-O,  122,  123 

214,  215,  218 

IV  Measures  of  performance:   MES103,  104,  109 

V  Results: 

1.  Next  to  shortest  imminent  operation  rule  (RUL105)  least  work 
remaining  rule  (RUL122)  was  the  best  as  far  as  percentage 
ULilization  of  shop  capacity  is  concerned. 

2.  Shop  capacity  utilization  is  minimum  for  most  work  remaining 
rule  (RUL123)  at  medium  load. 

3.  Rules  which  depend  on  monetary  value  of  jobs  (RUL214j  215) 
are  the  poorest  at  light  and  medium,  load  as  far  as  percentage 
utilizati.on  of  shop  capacity  is  concerned.   However,  they  v;ere 
comparable  at  heavy  loads. 

A.   It  is  expected  tliat  three-class  monetary  rule  (RUL216)  v.'ould 

be  better  than  two-class  monetary  rule. 
5.   Mean  lateness  of  jobs  (MES103)  is  decreased  as  the  work  content 

of  the  subsequent  queue  increases  before  the  queue  becomes 

critical  (RUL3.11)  .   That  is,  the  mean  lateness  decreases  as 

parameter  a  increases. 
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EXPERIMENT  3 
INVESTIGATOR         Kurtani  and  Nelson  [40] 
COMPUTER/PROG.  Li\NG  IBM  709/-- 

I   bilOP  CMARACTERISTICS 

1.  Shop  structure 

2.  Shop  type 

3.  LN'Unbcr  of  groups 

4.  Nurrber  of  machines /group 

5.  Shop  utilizatioa 
II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals 

2,  Distr.  of  job  routings 


Net  work  of  queues 
Job  shop 
4  and  8 

1  in  each  group 


Polsson  and  Erlang  (k  =  2) 


Known  probabilistic  transition 
matrix 


3.  DLstr,    of   Number  of  - 
operatioiis/job. 

4,  Distr^of   processing    time        Exponential   and  Erlangs    (k  -   2) 


256   runt 
Empty 


5.  Number  of  runs 

6.  Initial  Shop  State 

7.  Initial  number  of 
released  jobs 

8  .   Ru  n  -  i  a  1  e  a.e  t  Ii 


9.      Number   of   jobs    considered      1000 
for   collecting   statistics 

10,  J  Ui  n-  o  u  t   1  eng  th 

11,  Jobs/run  or  simulation     1000  jobs 
time 

12,  Due-date  multiplier        1,8  and  17 
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EXPF.RIMt:NT   3    (cent.) 
Ill      Scheduling   rules:      RUL102,    105 
IV     Measures   of   perfonnonce:      MES103,    109 
V     Results: 

1.  First    in   queue   first   served   rule    (RUL103)    riiinimizes    the 
inaxi.irum   flow    time   of   job, 

2,  Shortest   iivmiinent   operation  rule    (RUL105)    mininizes   average 
flow   time   of   jobs. 
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EXPERIMENT  4 
INVESTIGATOR  Conway  and  Ma.well  [IS] 

COMPUTER/ PROG.  L.\MG. 

I   SHOP  CHARACTERISTICS 

1.  Shop  structure 

2.  Shop  type 

3.  Number  of  groups 

4.  Number  of  machines /group 

5.  Shop  utilization 
II   EXPERIMENTAL  CONDITIONS 

1.   Distr.  of  job  arrivals 


Network  of  single  queues 
Job  shop  and  flow  shop 
2 ,  3  ,  6 ,  and  9 
1 


2.  Distr.  of  job  routings 

3.  Distr.  of  Number  of 
operations /job , 


Number  of  jobs  in  shop  were  kept 
constant  =2,  4,  6  times  number 
of  machines 

Uniform 

Geometric,  mean  =  M,  9  maximum 
operations  (Job  shop) 


4.  DLstr.  of  processing  time   Exponential,  mean  =  10  time  units 

5.  Number  of  runs 

6.  Initial  Shop  State  Empty 

7  .   I  111  t  i  al  numb  e  r  o  f  - 

released  jobs 

8.  Run-in  length  - 

9.  Number  of  jobs  considered   2000 

for  collecting  statistics  , 

10,  Run~o'.'.t  length  - 

11.  Jobs/run  or  simulation      2000  jobs 
time 


12.   Due-date  multiplier 
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EXPERIMENT  4  (cout:.) 
Ill   SchaJuliag  rules:   RUL106,  108,  111,  221 
IV  Meajarcs  cl  perf  orrn.ance :   ML.SIOS,  109,  110,  213,  216 
V  Results: 

1.  There  is  little  difference  between  perf ormance  in  pure  job 
shop  and  that  in  pure  flow  shop, 

2.  The  transition  matrix,  v;hich  determines  the  probability  of 
subsequent  operation  of  a  job  on  a  specific  machine,  is  not 
a  principal  determinant  of  priority  rule  performance. 

3.  Longest  imminent  operation  rule  (RUL109)  maximizes  what  the 
shortest  imminent  operation  rule  (RUL105)  minimizes;  however, 
both  the  rules  have  high  variance  of  flow  time. 

4.  Subsequent  operation  rule  (RULlll)  exhibited  considerably 
better  performance  than  the  random  rule  but  it  did  not 
surpass  the  shortest  imminent  operation  rule. 

5.  Advantage  secured  in  mean  shop  flov;  tim.e  of  job  is  directly 
proportional  to  the  fraction  of  time  the  shortest-operation 
time  rule  is  used  in  the  compound  rule  of  alternating  shortest 
operation  and  first  come  first  served  rule  (RUL221) . 

6.  In  a  system  consisting  of  a  network  of  queues,  there  is  a 
considerable  experimental  support  to  the  conjectures  that: 

6.1  All  local  priority  rules  such  that  each  priority  class 
has  some  expected  processing  time,  are  equivalent. 

5.2  The  shortest  operation  lule  (RUL105)  is  optimal  with 
..■•i-p^.ct  to  the  set  of  all  local  priorit/  rules. 

7.  The  shortest  operation  rule  (.RuLlOf-)  appeared  to  be  highly 
insensitive  to  errors  of  estimating  processing  times. 
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EXPERIMENT  5 
INVESTIGATOR  Gere  [29] 

COMPUTER/PROG.  LANG.     IBM  709/Foi-'i:rau 
I   SHOP  CHARACTERISTICS 

1.  Shop  structure  Network  of  single  server  queues 

2.  Shop  type  Job  shop 

3.  Number  of  groups  4  to  16 

4.  Nuraber  of  machines /group      1 

5.  Shop  utilization  - 
II   EXFERIMENTilL  CONDITIONS 

1.  Distr.  of  job  arrivals        - 

2.  Distr.  of  job  routings        - 

3.  Distr.  of  Number  of  - 
operations /j  ob . 

4.  Distr.  of  processing  time  Rectangular,  1  to  10  hrs.  per  operation 

5.  Number  of  ruias 

6.  Initilal  Shop  State         Empty 

7.  Initial  number  of 
released  jobs 

8.  Run--in  length  - 

9.  NuAiber   of   jobs   considered      1.      25   files   each  having    6-20   jobs, 
for   collecting   stiitistics      2.      16   files   each   having    20-60    jobs. 

10.  Run-out   length 

11.  Jobs/run  or  simulation     1.   Static  case:   25  files  each  with 
tine  6-20  jobs,  1-16  operations 

2.   Dynaiaic  case:   16  files  each  i;vith 
20-60  jobs,  1-16  operations. 

12.  Due-date  muJ  tiolier        2  to  4 
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E>:PERIMENT    5    (cont.) 

Ill      Scheduling   rules:      RULlOl,    103,    105,    llA,    207,    208,    209--212,    226, 

228,    229 

IV  Measures  of  performance:   MES102,  .105,  115,  116 

V  Results: 

1.  Priority  rules  which  considered  job  slack  time  (RUL114,  226, 
228  and  229)  were  more  effective  than  the  rules  in  v/hich 
priority  is  given  at  random  (RULlOl  and  103) , 

2.  If  all  the  jobs  have  different  number  of  operations, 
remaining  job  slack  ti.ne  rule  (RUL114)  was  somewhat  better 
than  remaining  job  slack  per  operation  rule  (RUL226) . 
However,  \vhen  each  job  has  sa:ne  number  of  operations  there 
is  no  significant  difference  between  the  two  rules. 

3.  Ratio  of  modified  job  slack  tine  to  tiuie  uui:il  due  date  rule 
(RUL229)  was  no  more  effective  than  RULIl't,  for  different 
machine  loading. 

A.   For  both  static  and  dynaraic  problems  job  slack  time  rule 
(RULllA)  was  significantly  more  effective  t\\a:n   shortest 
imminent  operation  rule  (R(jr-.105)  . 

5.  Non-random  rules  are  signif icaatly  more  effective  than 
random,  rules. 

6.  There  is  little  difference  in  ef f ectiveriess  between  the 
rules  v.hich  use  som.e  property  of  job  slack  time. 

7.  The  alternate  operation  and  look  ahead  heuristics  are  effective, 
both  Individually  and  collectively. 
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EXPEPTMENT  6 
INVESTIGATOR  LeGraude  ['O] 

COMI'UTER/PROG.  LANG.       IBM  7090/- 
I   SHOP  CHARACTERISTICS 

1.  Shop  Structure  Labour  restricted  fabrication  shop. 

2.  Shop  type  Job  shop 

3.  Number  of  groups  115  machine  groups  &  47  labour  classes. 

4.  Number  of  iriachines /group   Total  1000  machines  &.   400  to  500  raen. 

5.  Shop  utilization  Labour  utilization  of  60-70% 
II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals     KnoxvTi  historic  data 

2.  Distr.  of  job  routings  - 

3.  Distr.  of  Number  of  - 
operations /job . 

4.  Distr.  of  processing  time  Negative  exponential,  mean  -  mean 

processing  time  of  machine  group 
concerned 

5.  Nuraber    of    runs  - 

6.  Initial  Shop   Sl:ate  - 

7.  liiitial  number   of  - 
released  jobs 

3,       I'vua-in   length  - 

9.      Number   of   jobs   coiisidered     3000  jobs 
for   collecting   statistics 

10.  R'.nv-out    length 

11.  Jobs /run  or  simulation     3000  jobs 

('.ime 

12.  Due-date  multiplier 
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EXPERIMENT  6  (coul:.) 

Ill   Scheduling  rules:   RULlOi,  103,  105,  115,  116,  226 

IV  Measures  of  performance:   MESllO, 

202,  214,  217 

V  Results: 

1.   The  order  in  v.'hich  the  rules  were  effective,  starting  with 

the  'best'  rule,  was  as  follows: 

1.  Shortest  iraniinent  operation  (RULlOo) 

2.  Ratio  of  reaiaiuing  job  slack  to  total  remaining 
number  of  operation  (RUL226) 

3.  First  in  queue  first  served  (RUL103) 

4.  Earliest  start  date  (RUL115) 

5.  Earliest  job  due  date  (RUL116) 

6.  Random  rule  (RULlOl) 

2.  The  above  rating  was  made  under  the  v;eighting  systems  which 
assigned  weights  to  the  different  measures  of  effectiveness. 
The  job  completion  meiisures  were  given  the  highest  weighting 
factors. 
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EXPKRIMENT  7 
INVESTIGATOR  Conway  [12,16] 

COMPUTER/rROG .  l.ANG  .       IBM  7090/Siiuscript 
I   SHOP  CIIARACTERISTiCS 

1.  Shop  structure  Network  of  single-server  queues 

2.  Shop  type  Job  shop 

3.  Number  of  groups  9 

4.  Number  of  iuachines /group  1 

5.  Shop  utilii^ation  90%  . 
II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals  Poisson 

2.  Distr,  of  job  routings  Uniform  (1,M) 

3.  Distr.  of  Number  of  Geometric,  p  -   1/M 
operations /job  . 

4.  Distr.  of  processing  time  Exponential,  y  -  1 

5.  Shop  utilization  130 

6.  Initial  Shop  State  Empty  shop 

7.  Initial  number  of  50 
released  jobs 

8.  Run-in  length  400 

9.  Number  of  jobs  considered  8700  (401  to  9100  serially) 
for  collecting  statistics 

10,  Run-out  length  900 

11,  Jobs/run  or  simulation  10,000  jobs 
time 

12,  Due-date  multiplier 
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EXPERIMENT  7  (cont:.) 

in   Scheduling  rules:   RULlOl  -  106,  108  -  130, 

20/'r,  214  -  219,  221  --  226,  231  -  236 

IV     Moasures   of   pari omance:      RES201,    204,    207,    210 

V     Rosultis: 

1.  The  shortest  imminent  operation  rule  (RUL105)  clearly 
dominates  all  the  other  rules.   Combination  of  shortest 
operation  rule  with  any  other  rules  seems  always  beneficial. 

2.  Shortest  imminent  operation  rule  is  not  the  best  rule  with 
any  of  the  work  content  measures  (ri£S204  to  MES212)  . 

3.  The  v7ork  content  measures  rules  involving  processing  time 
factor  are  the  best. 

4.  For  three  out  of  the  four  methods  considered  for  settiv.ig  due 
dates,  shortest  imminent  operation  rule  caused  fewer  jobs 

to  be  late  than  any  other  procedure. 

5.  The  variance  of  distribution  of  shop  times  is  smaller  for 
shortest  imminent  operation  rule  than  for  aa^'  other  basLc  i-ule, 
except  first  in  shop  first  served  rule  (RUL102) . 

6.  Shortest  imminent  operation  rule  is  not  sensitive  to  errors 
in  estim.ating  processing  t  Lm.es. 

7.  Rule  in  which  priority  is  assigned  according  to  the  ratio  of 
remainii'ig  slack  to  remaining  number  of  operations  (RUL226) 
V7as  clearly  super  i.or  to  the  other  due  date  dependent  rules. 

8.  Where  due  dates  are  assigned  at  the  time  of  job  arrival, 

an  allowance  p.copcctional  to  the  total  processing  time  appeared 
to  be  better  than  the  other  methods. 
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EXPERIMENT   8 
rWVESTIGATOR  Carrol    [11] 

COMPUTE!? /FROG.    LANG. 

I      SHOP    CliARACTERISTICS 

1.  Shop   structure  Assembly   shop 

2.  Shop    type  Job   shiop  having   single   and  raultiple 

component  jobs 

3.  Number  of  groups  8 
A.  Nuniber  of  machines/group  1 
5.      Shop  utilization  80% 

II      EXPERIMENTAL   CONDITIONS 

1.  Olstr.    of   job   arrivals  Poisson 

2.  Dlstr.    of   job   routings  - 

3.  Dlstr.    of  Number   of 
operations / j  ob , 

4.  Dlstr.    of    processing    time     Exponential 

5.  Number  of   runs 

6.  Initial  Shop   State  - 

7.  Initial  number   of  - 
released   jobs 

8.  Run-in   length  - 

9.  Number   of   jobs    considered  - 
for   collecting   statistics 

10.  Run-out   leiv-th 

11.  Jobs/run  or   simulation  1.      Single   component  runs:      3000 
time  orders   involving   A0,000   tasks 

2.      Multiple   components   runs:      2000 
orders    involving   30,000   to   60,000 
tasks 

12.  D'lC-date  multiplier  - 
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EXPERIMENT  8  (coat.) 

HI   Sclieduling   rules:      RUL102,    103,    105,    108, 

226,    230 

IV  Measures  of  performance:   MES103 

V  Results: 

1.  Mean  jol)  tardiness  (MES106)  for  COVERT  rule  was  less  as 
compared  to  truncated  shortest  iromlaent  operation  rule 
(RUL108) 

2.  V'ith  COVERT  the  distribution  of  job  lateness  is  skewed  such 
that  jobs  completed  on  the  due  dates  are  the  mode  of  the 
distribution  with  very  few  late  jobs. 

3.  At  different  shop  utilisation  levels  and  with  dLffereat 
flow  allowance  there  was  no  significant  difference  in  the 
performance  of  COVERT  rule. 
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FXPKRTXFNT  9 
INVESTIGATOR  Orkin  [52] 

COMPUTER/PROG.  LMG.   CDC  -  .16C4/Fortraa  63,  Codap  -  1 
I   SHOP  CHARACTERISTICS 

1.  Shop  siiructure  Network  of  single  server  queues 

2.  Shop  type  Job  shop 

3.  Number  of  groups  8 

4.  Number  of  rfiachines/group    1 

5.  Shop  utilization  'J0% 
II   EXPERIMENTAL  CONDITEONS 

1.  Distr.  of  job  arrivals      Interarrival:   Geonetric 

2.  Distr.  of  job  routings      Uniform 

3.  Distr.  of  Number  of         Geometric,  p  -  1/8 
operations/job . 

4.  Distr.  of  processing        Geometric,  p  "  i/4 
time 

5.  Number  of  runs  - 

6.  Initial  Shop  State  Empty 

7.  Initial  number  of  15 
released  jobs 

8.  Rup.-in  length  215 

9.  Number  of  jobs  considered    6000 
for  colJ acting  statistics 

10.  Run-out  length  - 

11.  Jobs/run  or  simulation      6215  jobs 
time 

12.  Due- date  multiplier         5  and  9 
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EXPERIMENT  9  (cent.) 

Ill   Scheduling  ruJes;   RUL103,  105,  113,  116,  117, 

203,  226,  239 

IV     Measures   of   performance:      MES10  5,    106,    107,    109,    1.10 

V      Rcijults: 

1.   Setting  operation  clue  dates  proportional  to  processing  time 

results  significantly  lover  values  of  mean  job  tardiness 

than  setting  equally  spaced  operation  due  dates. 
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EXPERIMENT   iO 
IisiVESTIGATOR  Oldziey    [31] 

COI€'UTER/PROG .    LANG.      CDC  -    IGOA/T.PE,    Codap   1 
I      SliOP    CHARACTERISTICS 
1.      Shop   structure 


2.  Shop  type 

3.  Nunbar  of  groups 

4.  Number  of  machiaes/group 

5.  90% 

II   EXPERIMENTAL  CONDITIONS 

1.  DisLr.  of  job  arrivals 

2.  Dlstr.  of  job  routings 

3.  Di-str.  of  Number  of 
operations / j  ob . 

4.  Distr.  of  processing 
time 

5.  Number  of  runs 

6.  Initial  Shop  State 

7.  Initial  number  of 
released  jobs 

8 .  Ru  n  -  L  n  1  e  n  ■  <  t  h 


Network  of  single  server  queueing 
process 

Job  shop 

8 

1 


Geometric  inter arrival 

Uniform 

Geometric,  p  =  1/8;  max  of  63 
operations  per  job. 


Geometric, 


1/4 


Empty 

15  (for  sample  size  of  6v000)  , 
50  (for  sample  size  of  2000) 

200 


9.   Number  of  jobs  considered  2000  (201^^  to  2200"'^)  and  6000 

for  collecting  statistics  (20ist  to  6200'-^"^) 

10.  Run-cut   length  300 

11.  Jobs/run   or   simulation  2500  jobs    (fov/  factor  =   6),    6500 
time  (flov/   factor  -    5) 


12.      Due-date  multiplier 


4,5,6 


EXPERIMENT  10  (coni:.) 

Ill   Scheduling  rules:   RULlOl,  103,  105,  116, 

201,  202,  203,  206,  237,  238 

IV     Measures   of   perfomauce:      IIES103,    lOA,     LG6,    107,    110 

V     Results: 

1.  Sun  of  operation  slack  time  and  weighted  operation  process iiig 
time  rule  (R1JL238)  with  weighting  factor  of  15  reduced 

mean  job  tardiness  to  half  that- of  shortest  operation  rule 
(RULi05)  v/hile  reducing  variance  of  lateness  of  jobs  to 
less  than  1/3  the  value  obtained  v/ith  shortest  operation 
rule. 

2.  RUL201  improved  the  measure  of  mean  tardiness  marginally 
over  that  of  RUL238. 

3.  RUL239  showed  better  performance  over  RUL201  and  RUL23S  by 
the  inclusion  of  load  smoothing  factor?  cO .  T  t ( t) . 
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EXPliRIMENT   11 
INVESTIGATOR  Ellon   and   Hodgsoa    [23] 

COMPUTER/PROG.    LANG. 

I      SHOP    ClIAP^VCTERISTICS 
1.      Shop   structure 


2.  Shop  type 

3.  Number  of  groups 

4.  Number  of  machines/group 
3.   Shop  utilization 

II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals 

2.  Distr.  of  job  routings 

3.  Distr.  of  Number  of 
opera.':  ions /job . 

4.  Distr.  of  processing 
time 

5.  Number  of  runs 

6.  Initial  Shop  State 

7.  Initial  number  of 
released  jobs 

8 .  Ru n- i n  1 c  n e  th 


A  queing  system  of  two  identical 
servers  in  parallel 

Parallal  processor  shop  with  single 
queue 

Tv.'o  in  parallal 

1 

p  =  0.5,  0.6,  0.7,  0.8,  0.9 


Poisson 


Each  job  requires  one  operation 
Negative  exponential  distribution 


Emp  ty 


100 


9.      N'.imber  of  jobs   considered      1800 
for   collecting  statistics 


10.      Run-out   length 


100 


11.  vjobs/run  or  simulation  2000  jobs 
time 

12,  Due-date  multiplier 
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EXPERIMENT  11  (cont.) 

Ill   Scheduling  rules:   RULlOl,  103,  105,  109,  116 

IV  Measures  of  perf onnance:   MESlOl,  107,  108,  111,  113,  114,  117, 

201,  216 

V  Results: 

1.  Willi  respect  to  the  measure  of  mean  and  standard  deviation  of 
job  waiting  times,  it  v/as  shov;n  that  shortest  iroiriineut 
operation  rule  perfornied  the  best  and  longest  imminent 
operation  rule  (RUL109)  perforTned  the  poorest. 

2.  For  job  due  date  multiplier  of  8,  the  results  of  job  due 
date  rule  (RUL116)  was  close  to  shortest  operation  rule. 

3.  The  performance  measures  of  job  flow   time,  facility  idle  time 
and  job  delay  factor  are  independent  of  loading  rules  but 
dependent  on  the  load  ratio  (p )  of  the  system. 

4.  Shortest  operation  rule  was  considered  the  best  for 
minimizing  measures  of  job  waiting  times,  flow  times,  machine 
idle  time,  delay  factor  and  queue  lengths. 
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EXPERIMENT  12 
INVESTIGATOR  Moodio  nud  Roberl-.^  [AS] 

COMPUTER/PROC.  LANG.   IBM  7094/GPSS  III 
I   SHOP  CHAP^ACTlRISTICS 

1.  Sb.op    structure 

2.  Siiop    type 

3.  Number  of  groups 

4.  Number  of  machines /group    1 

5.  Shop  utilization  - 
II   EXPERIMENTz\L  CONDITIONS 

1.   Distr.  of  job  arrivals 


Identical  s e rv e v s  wi th  single  q u eu e 

Parallal  processor  shop 

24 


2.  Distr.  of  job  routings 

3.  Distr.  of  Number  of 
operations/job . 


Known  historic  distribution, 
arrival  rate  =  10,  11,  12 


Number  of  operations  per  job  is 
cons is  tan t 


4.  Distr.  of  processing  tine   Known  historic  distribution 

5.  Number  of  runs 


6.  Initial  Shop  State 

7.  Initial  number  of 
released  jobs 

8.  Run-in  length 

9.  Number  of  jobs  considered 
for  collecting  statistics 

10.  Run-out  length 

11.  Jobs/run  or  simulation 
time 

12.  Due-date  multiplier 


Emp  ty 


4S0   orinutes   per  day;    365   days 
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EXPERIMENT  12  (couC.) 

Ill   Scheduling  rules:   RUL103,  lO^,  114,  116,  240 
IV  Measures  of  perfo-niiance:   MES102,  103,  104,  106,  107 
V  Results: 

.1.   Job  due  date  rule  (RUL116)  failed  to  minjiuize  the  maximum 
lateness  in  a  parallal  processor  shop, 

2,  The  weighted  objective  rule  (RUL240)  showed  remarkable 
consistency  for  changing  shop  loads. 

3.  Mean  job  tardiness  seemed  to  be  best  minimized  by  remaining 
job  slack  time  rule  (RUI,114)  and  weighted  objective  rule, 
but  for  highly  loaded  shop  shortest  operation  (RUL105)  and 
weighted  objective  rule  were  the  best. 
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EXPERIMENT  13 
INVESTIGATOR  Neimeir  [50] 

COMPUTER/PROG.  L/\NG .   CDC  1604/CTP  (Cornell  list  prooessor) 
I   SHOP  CHARACTERISTICS 

1.  Shop  structure 

2.  Shop  type 

3.  Number  of  groups 
A.   Number  of  machines /group 
5.   Shop  utilization 

II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals 

2.  Distr,  of  job  routings 


Netv/ork  of  single  queues 

Job  shop  with  alternating  routing 

9 

1 

90% 


3.  Distr.  of  Nuaiber  of 
operations/ job . 

4.  Distr,  of  processing 
time 

5.  Number  of  runs 

6.  Initial  Shop  State 

7.  Initial  number  of 
released  jobs 

8.  Run-in  length 


Poisson,  mean  interarrival  =  1.09 

Uniform 

Geometric,  p  -  1/9 


Exponential,  mean  -  0.96 


Emp  ty 
50  jobs 


40,  200,  400  time  units  corresponding 
to  short,  medium  and  long  run 


9.   Number  of  jobs  considered  After  every  30  t  v.:\q   units 
for  collecting  statistics 


10.   Pom-out  length 


11,   Jobs/run  or  simuia;;lon     1200,  5000,  10000  time  units  cov- 
tirac  responding  to  short,  medium  and 

long  run 


12.   Due-data  multiplier 
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EXPERIMENT  13  (corit.) 

Ill  Scheduling  rules:   RiILl02,  105 

IV"  Measures  of  perf on^uuice:   MES108, 

201,  216 

V  Results: 

1.  Alcernative  routing  could  be  prot'italjlj'  employed  in  a  job 
shop  under  first  in  queue  first  served  rule. 

2.  There  was  greater  improvement  in  performance  under  the  first 
coi:;e  first  ser^/ed  rule  (RUL103)  than  under  shortest  operation 
rule  (RUL105) J  with  introduction  of  alternative  routing. 

3.  Alternative  routing  vjith  the  selection  of  machine  with  least 
number  of  vjaiting  jobs  had  slightly  better  performance  than 
that  of  selecting  the  machine  with  least  amount  of  work  in 
queue. 
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hXi'LlUiMENT  14 
INVESTIGATOR  Thompson  [CO] 

COMPUTER/ PROG.  LANG.      CHC  3 400 /GASP  iL 
I   SHOP  CHAl<i\CTERISTICS 

1.  Shop  st;ructure  Assembly  job  sliop 

2.  SViup  type  Job  shop 

3.  Number  of  groups  20 

4.  Number  of  machines/group    1 

5.  Shop  util.i/.ation  75% 
II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals     Poisson  A  =  20  products  per  day 

2.  Distr.  of  job  routings 

3.  Distr.  of  Number  of        Triangular,  mean  -  15,  range  =  10 
operations/job .  to  20 

4.  Distr.  of  processing  time   Exponential,  m.^an  ~  0.1  day 

5.  Number  of  runs  37 

6.  Initial  Shop  State         Empty 

7.  Initial  number  of  - 
releaocd  jobs 

8.  Run-in  length  100 

9.  Number  of  jobs  considered  300 
for  collecting  statistics 

10.  P.un-out  length 

11.  Jobs/run  or  simulation     400  jobs 
time 

12.  Cue--date  multiplier 
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EXPERIMENT  14  (coat.) 

Ill   Scheduling  rules:   RUL102,  105,  119,  206 
IV  >'e;isures  oJ:  performance:   MES108 
V  Results: 

1.   Performance  of  maximum  number  of  operation  rer.aining  rule 
(RUL119)  v;as  as  v;ell  or  better  than  any  of  the  rules  tested 
for  assembly  job  shop,  however  no  cptimality  can  be  guaranted, 
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EXPERIMENT    15 
INVESTIGATOR  EiJ.on  and   Coi-.r.erill    [22] 

COMPUTER/ PROG.    LANG.  IBM   7090/- 

I      SHOP    CHAPvACTERISTICS 
L.      v'ihop   structure 

2.  Shop,  type 

3.  >!umber  of  groups 

4.  Number  of  machines /group 

5.  Shop  utilization 
II   EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  iirrlvals 

2.  Distr.  of  job  routings 


Network  of  queues 

Job  shop 

4 

1 

p  =■-■  0.4,  0.5,  0.6,  0.7,  0.3 

Poisson 
Random 


3.   Distr.  of  Number  of 
cpcrations/j  ob . 


Rectangular  distribution  number  of 
operation 


4.  Distr.  of  processing  time  Negative  exponential   per  job  tv/o  to 

four 

5.  Number  of  runs  - 

6.  ItiitJal  Shop  State         Empty 

7.  Initial  number  of  - 
released  jobs 

8.  Pom- in  length  - 

9.  Number  of  jobs  considered   5000  jobs 
for  collecting  statistics 

10,  Run-out  length  - 

11.  Jobs/run  or  c.  ^'v.ul  ■  r.ion     5000  completed  jobs 
time 


12.   Due-ditc  multiplier 


47 


EXPERIMENT  15  (cont:,) 

III      Scheduling   rules:      RUL103,    105,    107-109.    116,    122-124 

206,    210,    220 

IV  Measures  o£  perfonncrnce:   KESlOl,  108,  ill,  113,  114 

V  Results: 

1.  The  difference  between  performance  of  various  loading  rules 
b<?.come  increasingly  pronounced  as  the  load  ratio  increases, 
the  shortest  operation  rule  (RUL105)  being  more  effective 
than  first  come  first  served  (RUL103)  and  longest  operation 
rule  (RUL107) . 

2.  The  shortest  operation  rule  performs  best  with  several 
criteria,  but  not  for  minimizing  the  variance  of  flow  times 
or  lateness  of  jobs. 

3.  The  total  simulation  run  tine  (for  completing  50C0  jobs) 
and  machine  idle  time  depend  on  the  load  ratio  (p)  but  not 
on  the  loading  rules. 
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CHAPTER  III 


SI^rULATION  MODEL 


This  chapter  deals  v;ith  the  discussion  of  the  job-shop  model  on 
which  our  experiments  are  based.   The  system  description,  experimental 
cop.dLtions,  scheduling  rules  and  the  measures  of  perf omumce  are 
described  in  detail  below.   The  simulation  computer  nrogr^im  together 
with  the  purpose  of  the  subroutines  and  the  definition  of  the  program 
variables  is  presented  in  Appendix  C.   A  general ,  flow  chart  of  the 
job-shop  simulation  is  given  in  Figure  3.1. 

3.1  Model  Cescription 

The  job-shop  ccnsidered  in  this  paper  is  a  network  of  single-server 
queues  in  which  a  job  may  join  all  or  a  subset  of  queues,  waiting  to 
be  served  at  different  machine  groups.   Each  machine  group  is  an  in- 
depei;;rier;t  work  center.   All  machines  in  a  particular  group  are  identical 
as  fa'i"  as  V.ie    type  of  v.-ork  to  be  performed  on  a  job  is  concerned.   In 
other  words,  jobs  x\'aiting  at  any  machine  group  form  a  single  queue  and 
are  processed  on  any  of  the  m£ichine  which  becomes  idle  first. 

In  operating  such  system,  the  utilization  of  the  shop  may  be 
assumed,   Tlic  utilization  of  the  s'lop  may  be  discussed  as  follows.   Each 
Lime  a  job  is  processed  on  a  machine,  its  services  are  utilized  and 
the  machine  is  said  to  be  busy,  otherwise  it  remains  :Ldle.   Hence,  the 
utilizat:ion  of  a  machine  can  be  defined  as  thie  ratio  of  the  total  busy 
time  of  that  machine  to  th>e  total  time  for  \/hich  it  has  been  available. 
TVve  uti].i/:ation  of  thi^.  s^iop,  as  a  whole,  is  the  average  oi    the,   utilizationf 
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1 

3 
4 
5 
5 


St; 


"■'3 


Zi 

Input: 
Number  of  raachina  groups 
Number  of  machines  per  group 
Job  sample  size  or  simulation  time 
Number  of  simulation  runs 
Scheduling  rule 
Due-date  procedure 


Job  File: 

1.  Inter-arrival  times 

2,  Number  of  operations 

3,  Machine  ordering  (routing) 

4.  Processing  times 


1 


Shop  Simulation: 

1.  Pvclease  the  job  to  the  shop 

2.  Schedule  the  job  to  available 

machine  according  to  its 
routing 

3.  Collect  statistics  for  the 

individual  job 

4.  Release  the  job  out  of  the  shop 


No 


Job  s.-;!/:iple 
completed? 


Figure  3,1  Cenec'sl  simulatloa  flow  chart 
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of  individual  r:ac'.iines.   For  a  network  sysLcu  of  job -shop,  the 
utilization,  called  load  ratio  in  [22],  can  be  defined  such  that 

m 

where 

U   utilization  of  the  shop 

n   job  sample  size 

L   mean  processing  time  per  job 

M   total  number  of  machines  in  the  shop 

T   total  simulated  run  time 

In  summary,  the  shop  characteristics  considered  in  this  investi- 
gation are: 

1.  Shop  structure:  Netv/ork  of  single-queues 

2.  Shop  type:  Job-shop 

3.  Number  of  machine  groups:         9 
A.   Munber  of  machine  per  group:      1 

J.   Ctilization  of  shop;  ~  93% 

3.2  Experimental  Conditions 

In  order  to  perform  a  simulation  experiment,  certaxi\  coiiditions 
have  to  be  specified.   The  following  conditions  are  considered  In  our 
experiments. 

1.   Job  Arrivals ;   At  a  time,  only  one  job  is  released  to  the  shop. 
The  inter-arrivals  of  jobs  are  obtained  from  an  exponential  distribution, 
However,  Erlang  distributions  vj  itli  different  pai'ameters  are  also  used  in 
other  experiments.   The  exponential  m.ean  inter-arrival  time  was  fixed 
at  1.0  time  units,  similar  to  [12]  so  that  coinparisoa  of  the  results 
can  '^e  "nt^o. 
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2,   Job  Routing:   On  arrivals,  each  job  is  equally  likely  to  start 
its  first  operation  on  any  of  the  raacliines  ~.v.   the  shop.   Furth.err.iore, 
if  the  job  has  more  than  one  operation,  it  ba??  equal  probability  of 
starting  its  next  operation  on  any  of  the  remaining  TP.achines  in  th.e 
shop.   Therefore,  the  machine  ordering  of  each  job  is  generated  randomly 
from  a  uniform  distribution  [1,9].   In  other  \Jords,  all  elements  of  the 
machine  ordering  matrix  have  equal  probability  value  1/9. 

■^ '      Number  of  Operations  Per  Job;   Since  each  job  may  be  processed 
more  than  once  on  a  certain  machine,  the  number  of  operations  for  each 
job  has  to  be  determined.   In  our  experiments,  the  distribution  of  the 
number  of  operations  per  job  j  is  governed  by  the  follov/ing  distribution: 

K.-l 
P(KJ  =  p(l-p)  ^    ,   K^  =1,  2,  ... 

The  above  probability  function  represents  a  geometric  distribution  'vith 

mean  equal  to  1/p  or  9,  where  9  is  the  total  number  of  machines  in  the 

shop.   Each  job  is  allowed  to  have  at  the  most  12  operations. 

^*'      Processing  Times;   In  some  experiments  the  processing  times  of 
operations  are  obtained  from  an  exponential  distribution.   The  mean  of 
the  distribution  is  such  that  the  desired  shop  utilization  is  obtained. 
Ilovvever,  Erlang  distributions  with  different  parameters  are  also  used  in 
other  experiments, 

5.   Number  of  Runs:   Each  rule  is  tested  for  one  simulation  run. 
To  gain  more  kao;;lodge  about  the  effectiveness  of  any  rule,  ic    is  pre- 
ferred to  have  more  runs  each  of  -//riich  has  different  sets  of  jobs.   The 
lack  of  computer  time  available  is  the  Tain  factor  for  considering  only 
one  run  per  rule. 
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6.  InJLlal  Shop  State;   Before  conductni\g  any  simulation  experi- 
ment, it  is  necessary  to  place  the  jhop  in  yo.ne  more  realistic  initial 
state.   In  inxtiating  an  experiment,  one  of  the  raethccls  is  to  assume 
an  empty  sliop.   However,  this  directly  affects  the  leogth  of  the  tun, 
which  should  be  long  enough  to  obtain  representative  results  of  the 
stead}'  state  conditions.   A  more  desired  procedure  v/ould  be  to  start 
the  shop  \ntl\   steady  state  conditions,  after  they  have  been  known 
from  some  initial  experimental  results.   For  the  sake  of  simplicity, 

the  shop  considered  here  starts  empty  as  has  been  assumed  in  [22,  23,  25]. 

7.  Initial  Number  of  Release  Jobs;   There  are  no   released  jobs 
\/nen  the  shop  is  started;  however,  the  loading  of  the  shop  could  be 
accelerated.   Eilon  and  Cotteril  [22]  has  reported  that  the  initial 
number  of  release  jobs  lias  very  little  effect  on  reaching  the  steady 
state  conditions. 

^*   Run-in  Length;   Since  the  shop  starts  in  an  empty  state,  a 
certain  time  period  is  necessary  before  it  reaches  the  steady  state. 
Any  statistics  collected  during  this  transient  period  are  likely  to 
be  not  representative.   Unfortunately,  there  is  no  practical 
ana].ytical  method  available  by  which  steady  state  condition  could  be 
determined  befoirehand.   Hence  one  has  to  depend  on  the  initial 
experiments  vyhich  liave  to  be  conducted  mainly  in  order  to  gain  an 
insight  in  the  stability  of  the  system.   For  determining  the  steady 
state  cou'litions  for  our  experiments,  the  absolute  differences  between 
mean  inter-arrivals  and  mean  inter-departures  are  computed  sequentially. 
The  steady  state  is  considered  to  have  started  when  such  difference 
is  nearly  vanished.   The  number  of  jobs  af~..er  which  the  steady  state 
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condition   can  be  starLed   is    deteruiined   for   each   experiment:  by   plotLing 
the  graph   of   inter-arrival   and   intor-depax'ture   times   of   the  jobs, 
see  Cho-pt.jr   j-V, 

9.      Number  of  Jobs   Considered   for  Collecting  Statistics;      As 
mentioned   earlier,    for   consistent  results   one  should   run   the  sinuilation 
experiment   for   long   runs.      A  long   run  neutrall?:e3    the   effects   of 
assuraiiig  wrong   initial  shop   conditions.      To   comprornize  betv;een   the 
con.puter   time   and   the   consistency   desired   in  results,    statistics   are 
collected   at   every   100  jobs   departed   from   the  shop.      The  final  statistics 
are  also      ollected   at   the   end   of    the   run, 

^^ '      R^m-Out  Length;      Run-out    length   is   necessary    to  remove   the 
effect   of  uncompleted  jobs,    specially    in   the  rules  v/hich   give  preference 
to    the   last   operation  of    the  jobs.      Performance  of   svich   rules    is   :_;reatly 
affev-^ted,    If   at   the  end  of    the   experiment   there  are  raany  jobs  waiting 
for    their   last   operations.      Since   the   ruJ.es   considered   do   not   fall   into 
the   above   category,    no   run-out   length   is    allcv;ed   in   this   investigation. 
Also    this    is    supported  by  previous    investigations    [22,    25]    which  have 
showed    that   assumptions   of  no   run-out  period   are  not  serious    to   affect 
i'.he   accuracy   of    the   analysis. 

■'-'-•      Jobs  Per  Run:      To    limit    the   computer   time  v/ithin  30  minutes 
for   ea-h   run,    the  sample  size  of   each   run   is    1300   departed  jobs. 

^^'      I^ "'^.^.:" P ^^.^ g^  Multiplier;      Before  a  job   is   released   to    the  sliop, 
a  due-date   is   determined.      I*'cur  methods    of   assigning   such   due-date  have 
been  proposed  by   Conway    [12],      These   are; 
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1.  A.lJuv-ip.co   propor  I,:-!  oTi;\]    to   the    Lolal  work   to  be  done   on 
the   job 

rr 
^-. 

.J 
D.    -   r.,   +  9      )        t., 

Although  Conway  has  not  nif;ntioned  the  reason  for 
determiiiing  the  value  of  the  constant  9,  we  believe 
that  it  is  the  total  number  of  machines  in  the  shop, 

2.  Allowance  proportional  to  the  number  of  operations  on 
the  j  ob . 


D.  =  r.,  H-  8.883  K. 
3  jl  J 


3 .   C o ns  t  an  t  all owan c e 


D.  =  v.,   +  78.7985 


4.   Random  allowance 


D.  =  r       -1-  157.597X. 


where 


D.  due   date   of   job   j 

r.-.  ready   time  of   job   j    for  operation  1 

K.  the   number  of   operations   on  job   j 

t ,.  the  processing   time  of   job   j    for   operation  k 

•  X, .      a  random  variable   for   job    i 
3 

It   has  been   found   in    [12]    that   assigning   due-dates   by    tlie   first 
method  is   superior   to   any  other  methods.      In   this    paper,    a  new  method 
of  assigning   due-date  is   proposed  which   is   expressed  such  that 


K. 


D.  -  r..  +  K.   y    t.. 


It  differs  from  the  first  ivieLhod  in  the  sense  that  the  allowance  is 
proportional  to  the  Duiiiber  of  operations  as  veil  as  the  total  work  to 
be  done  on  the  job.   In  method  1.,  the  allowance  is  computed  by 
multiplying  the  total  work  content  of  the  jobs  by  the  nuTbe.r  of  machine 
groups  regardless  to  the  n-.'  ^er  -J  operations  to  be  performed  on  the 
-job.   This  proposed  method  is  considered  to  be  more  logical,  since 
for  the  sa-a  total  work  content;  the  job  which  has  a  greater  number  c'. 
operations  is  more  likely  to  be  late  than  the  one  which  has  fewer 
operations.   Accordingly,  the  job  with  the  larger  minber  of  operations 
should  be  given  the  greater  allovjance. 
3.3   Scheduling  Rules: 

Since  one  of  the  purpose  of  this  investigation  is  to  test  the 
method  of  a-.signing  due-dates  as  described  earlier,  due-date  rules 
are  considtced.   The  following  due-date  rules  are  tested  together 
with  the  f irst-in-queue  iirsc-served  and  shortest  imminent  operation 
time  rules. 

1.  Shortest  imrainenc  operation  time  rula  (RULIOS) :   Select  from 
the  queue  of  waiting  jobs,  the  job  which  will  occupy  the 
machine  for  the  'i^inimum  period  of  time, 

R.(t)  =  t.p. 
3       3(1 

2.  First-m-queue  first-served  rule  (RUL103)  :   Select  the  job  which 
has  the  earliest  arrival  time  at  the  queue, 

R.(£)  -  r.^  . 
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3.      Earliest   job   clue-date    ruj.e    (RUL116)  :      3eiect   first    the   job 
which  has    the   earliest  job   due-data, 

x-.(.0    =   D,    . 

A.   Remining  job  slack  time  rule  (RUL,114)  :   Select  first  the  job 

which  has  the  minimum  free  time  left  after  all  its  remaiuing 

operation  times  are  deducted  from  its  due-date, 

K. 
R.(t)  -  D.  -  t  -   )    t.,  . 

IC— l(- 

5.   Remaining  job  slack  time  per  remaining  operations  (RUL226) : 
Select  first  the  job  which  has  minim.ura  free  tim.e  left  per 
operation, 

K. 
D.-t-  )    t., 

I,  (,)  = ^LA . 

•^        K.  -  fe+1 
J 

3.4  Measures  of  Performance: 

The  choice  of  measures  of  performance  depends  mainly  on  the  rules 
to  be  evaluated.   As  the  last  three  rules  considered  are  job  due-date 
rules,  the  m.easures  of  performance  adopted  should  indicate  the  deviation 
of  jobs  cc-jipletion  time  from  the  pre-planned  due-date.   Lateness  of  job, 
job  flow  time  and  number  of  jobs  which  were  behind  their  due-date  are 
examples  of  such  measures  v/hich  are  considered  in  our  experiments. 
Shop  characteristics  such  as  utilization  of  shop  and  average  number  of 
jobs  in  the  shop  are  also  computed.   To  summarize,  the  follovjing  measures 
of  perfonnance  are  considered: 


57 


1.  Job  lateness 

2.  Average  shop  flow  time  of  job 

3.  Number  of  jobs  late 

A.   Average  number  of  jobs  in  the  shop 

5.   Utilizacioa  of  the  shop 
The  next  chapter  is  devoted  to  the  discussion  of  the  results  and  their 
analyses. 
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CHAPTER  IV 


SIMULATION  RESULTS 


This  chapter  is  devoced  to  the  results  of  three  different  experi- 
ments totalling  14  simulation  runs.   A  representative  results  of  a 
particular  siiivalation  run  is  presented  in  Appendix  B  along  with  the 
graphs  v;hich  show  the  convergence  of  the  inter-arrivals  and  inter- 
dfcparturas  in  order  to  determine  the  steady  state  condition  fot  each 
run.   The  results  of  all  the  sii.iulation  runs  pertaining  to  each  'Experi- 
ment are  summarized  seperately  in  a  tabular  form. 

4 , 1  Experiment  I 

The  purpose  of  this  experiment  is  to  compare  the  results  of  this 
investigation  to  the  results  obtained  by  Conway  [12]  under  nearly  the 
same  set  of  experimental  conditions.   The  two  main  differences  are 
(1)  the  consideration  of  a  saiaple  size  of  1000  jobs  instead  of  lOjOOO 
jobs;  and  (/.)  the  setting  of  the  upper  limit  of  the  number  of  operations 
per  job  as  12  instead  of  39.   Experiment  I  is  necessary  in  order  to 
compare  the  results  of  experiment  II  in  which  the  nev;  procedure  of 
determining  the  job  duo-date  is  evaluated.   For  convenience,  the  results 
obtained  from  Conway  [12]    for  the  five  scheduling  rules  under  test  in 
tills  paper  are  displayed  in  Table  4.1. 

In  order  to  determine  the  number  of  departed  jobs  at  which  the 
shop  reaches  the  steady  state  conditions,  the  grand  averages  of  the 
inter-arrivals  and  inter-departures  of  the  jobs  are  plotted  in 
Figures  B-1  through  B-5  for  each  rule  tested.   It  is  found  that 
regard I'ss  of  tlie  rutr:,  thij  shop  reaches  the  steady  state  conditions 


Table   4.1 


Results   Obtained   from  Conway's   Experiraent 
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Scheduling 
rule 

FCFS(RUL103) 

SHOPN(R[JL105) 

SLACK(RUL114) 

DDATE(RUL116) 

S/OPK(RUL226) 


Lateness 

Shop 

flov7  times 

Number 
of   jobs 
late 

n> 

/erage 

Mean 

Std.    Dev. 

Mean 

Std,   Dev. 

jmber   of 
jbs   in 
.leues 

-4.49 

41.1 

74.43 

75.8 

3898 

53.87 

"44.9 

53,7 

34.02 

48.2 

437 

23.25 

-13.13 

20.8 

65.79 

80.8 

1916 

51.25 

-15.53 

20.8 

63.72 

82.4 

1544 

49.45 

-12.79 

15.0 

66.13 

73,8 

323 

51.65 

after  approximately  300  jobs  have  departed  from  the  shop.  Hence  the 

statistics  are  collected  after  the  departure  of  299   job  and  until  1300 

jobs  have  departed  from  the  shop.   Table  4.2  presents  the  results  of  the 
five  rules  tested  in  this  experiment. 

Table  4.2 

Results  of  Experiment  I 


Scheduling 
rule 


Lateness 


Mean    Std.  Dev. 


__S hop  flow  tim e_s_ 
Mean    Std.  Dev, 


Number   Average 
of  jobs  number  of 
late     jobs  in 


the  shop 

FCFS (RUL103 

-36.34 

25.82 

9.26 

5.72 

7 

11,30 

St{0PN(RULi05) 

"37.30 

25.22 

8.39 

6.17 

0 

10,35 

SLACK(RUL114) 

--35.38 

23,12 

3.12 

6.G9 

0 

J  0.06 

DDATE(RUL116) 

-35.34 

25.41 

8.01 

6.72 

0 

9,90 

S/OPN(RUL226) 

-34.79 

2';.  48 

8..  51 

5.98 

10 

10.46 

Note  that  the  last  coluiTtn  ia  both  tables  is  different. 

Significant  results  have  been  obtained  from  the  above  experinent. 
These  are  summarized  as  below: 

1.  The  shortest  Imminent  operation  rule  (RUL105)  is  no  doubt  the 
best  of  all  the  rules  tested  in  reducing  the  lateness  of  a  job.   The 
negative  value  of  lateness  indicates  that  the  job  is  completed  earlier 
than  the  due-date,  see  Table  4.2. 

2.  The  slack  time  rule  (RUL114)  and  the  earliest  job  due-date 
rule  (RUT. 116)  are  equally  good  in  reducing  the  lateness.   The  difference 
between  the  mean  lateness  is  negligible  to  conclude  the  superiority  of 
one  over  the  other,  see  Table  4.2. 

3.  In  comparing  the  results  of  Tables  4.1  and  4.2,  it  is  revealed 
that  the  numbers  of  jobs  late  under  all  rules  tested  are  much  larger  in 
Conway's  experiment  than  those  in  ours.   This  is  because  the  maximura 
number  of  operations  per  job  is  39  compared  to  ours  which  is  12.   Hence 
due  to  this  truncation  the  mean  number  of  operations  per  job  is 
approximately  4.2  instead  of  the  theoretical  mean  of  9. 

4.  The  mean  shop  flow  time  of  the  jobs  under  first- come  first- 
served  rule  (RUL103)  is  much  higher  than  that  of  the  shortest  imminent 
operation  rule  (RIIL105)  as  has  also  been  shown  by  Conway  [12].   However 
in  this  rcspecu  the  earliest  job  due-date  rule  (F-UL116)  surpassed  all 
the  rules  tested.   Furthermore,  the  moan  shop  flow  time  under  the  job 
slack  tinie  rule  (RUL114)  is  very  close  to  that  obtained  by  the  earliest 
job  due-date  rule. 

5.  Tl'ie  standard  deviation  of  shop  flew  tiiiies  is  not  proportional 
to  the  mean  of  the  shop  flow  time.   This  is  also   shoxvTi  by  Conway  [321. 
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6.  Even   though    tlie  shortest:   iimniuent   operation   ru.le  is    the  best 
in  minimizing   the  inean  job   lateness,    it   keeps    long  jobs  waiting    for 
much   longer   times    in   the  queues    than  most   of    the  other  rules    tested. 
This    is   evident  from   the   fact    that    the  maximum  shop   flov;  time   encountered 
under    this    rule   is    che   second  highest    (the   first   being   under    the   slack 
time  per  regaining  number   of  operations) .      The   reason   for   this    is    that 
the   shortest   operation   rule  selects    the   shortest   job    first   frc:n   the 
queue   of   v/aiting   jobs.      The  maximum    time    that   a  job  waits    in    the  shop 

is   shc-.m  below   for   each   rule    tested: 

Scheduling   rule  Maximum   time   in   shop 
FCFS(RUL103)  27.43 

SHOPN(Rl]L105)  36.67 

SLACK(RUL114)  35.51 

DDATE(RUL116)  35,38 

S/OPN(RUL226)  38.18 

7.  As  mentioned   earlier,    the   condition  of   the  shop   is   reflected 
by    the  measure  of   the  average   number  of  jobs   in   the  shop.      The   due- 
date   rule  is   the  best   of   all  rules    tested   in  having   few  jobs   in   the 
shop.      The   shortest   imminent   operation  rule   ranks    tlie   third;   liowever, 
the  difference   is    insignificant   to   conclude   that   it    is    inferior    to    the 
due-da;:e  rcule.      Iha  first-come   first-served  rule  has   on   the  average 
the  niaxiuium  number  of   jobs  waiting   in   the   shop.      This   is    in  agreement 
Xs'ith    th.e  results    detained  in    [12], 
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4 . '?.     Experiment  II 

Tlie  purpose  for  coiiducting  Experinent  II  is  to  test  the  new  pro- 
cedure of  assigning  due-dates  as  proposed  In  Chapter  III.   All  the 
experimental  conditions  of  this  experiment  are  exactly  the  same  as 
those  of  Experiment  1,  which  serves  as  the  basis  for  comparison  of 
the  results  obtained. 

As  before,  the  graphs  of  inter-arrivals  and  inter-departures 
are  plotted  seperately  for  each  rule  in  Figures  B-6  through  B--10j 
Appendix  B,   These  graphs  show  that  the  new  procedure  of  assigning 
due-dates  have  no  effect  on  the  sb.op  in  reaching  the  steady  state 
conditions.   The  shop  is  completely  free  frorr.  the  initial  transient 
effects  v;hen  approximately  300  jobs  have  departed  froin  the  shop. 
Table  4.3  displays  the  results  obtained  from  this  experiment. 

Table  4.3 

Results  of  Experiment  II 


Scheduling 
rule 


Lateness 


Shop   flow   times 


Mean    Std.  Dev.    Mean    Std.  Dev, 


Number   Average 
of  jobs   number  q£ 
late     jobs  in 


tiie  shop 

FCFS(RUL103) 

-35.33 

40.08 

8.22 

5.22 

185 

10,13 

SROPN(RUL105) 

-36,04 

39,64 

7.60 

5.49 

78 

9.47 

SLACK(RTJL114) 

-35.63 

38.23 

8.23 

7.13 

40 

10 .  19 

DDATE(RUL116) 

-35.51 

38.15 

3 ,  19 

7 ,  05 

44 

10,12 

S/0PN(RUL226) 

-35. C 2 

39 .  87 

8.53 

5.92 

178 

10.49 
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A  summary  of  conclusions  diawn  is  given  bela^?: 

1.  In  general,  the  mean  lateness  of  jobs  uirde?:  the  nev  procedure 
of  assigning  due-dates  has  very  little  effect,  an  increase  of  ]cs3  than 
2%  on  all  the  rules  tested. 

2.  In  cor-iparlng  the  results  in  Tables  4.2  and  4.3,  it  is  observed 
that  the  numbers  of  jobs  completed  late  are  larger  in  Experiment  II  than 
those  in  Experiment  I.   This  shows  that  the  new  procedure  proposed 
produces  the  condition  of  tight  job  due-dates  in  comparison  to  the 
procedure  in  which  the  allowance  is  computed  as  fixed  multiple  (equal 

to  number  of  machine)  of  work  content  of  job,  see  Experiment  1. 

3.  The  standard  deviation  of  job  lateness  is  increased  by  60% 
Eore  than  that  obtained  in  Experiment  I. 

A,   The  mean  shop  flow  times  for  the  non-due-date  rules  is  reduced 
by  10%  than  those  obtained  in  Experiment  I. 

5.  The  due-date  rules  are  not  significantly  affected  by  the  nev; 
procedure  of  assigning  due-dates.   The  mean  of  shop  flow  times  are  in- 
creased by  less  than  1%. 

6.  The  due-date  rules  are  not  much  sensitive  to  the  new  procedure 
as  far  as  standard  deviation  of  shop  flow  times  is  concerned.   The 
slack  per  remaining  number  of  operation  rule  (RUL226)  is  the  least 
sensitive  among  all  due-date  rules. 

7.  The  non-due-date  rules  gained  much  in  reducing  the  standard 
deviation  of  the  distribution  of  shop  flow  tim.es.   The  reduction  is 
approximately  10%  in  both  rules. 
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4 . 3  Experiment:  III 

The  purpose  o^  this  experlr.icnt  is  to  study  the  effect  of  different 
processing  tirae  distributions  on  the  perfonriance  of  the  iibove  rules. 
The  tv/o  different  distributions  selected  are  Erlang  4  a.nd  Erlang  8.   Due 
to  the  lack  of  the  computer  time,  it  was  necessary  to  collect  the  final 
statistics  after  600  jobs  had  departed  frora  the  shop.   The  reraaining 
experimental  conditions  are  Lhe  sane  as  those  of  Experiment  I.   The 
graphs  of  inter-arrivals  and  inter-departures  are  plotted  in  Figures 
B-11  through  B-14.   It  is  observed  that  under  the  shortest  operation 
laile,  tiie  shop  reaches  steady  state  condition  at  approximately  the  same 
stage  as  it  did  in  the  previous  tv;o  experiments.   However,  the  steady 
state  can  not  be  predicted  in  the  case  of  the  slack  per  remaining 
operation  rule.   Tl\is  is  mainly  because  of  small  job  sample.   Tlie 
results  of  this  experiment  are  suu-miari^^ed  belov/: 

Table  4.4 

Results  of  Experiment  III 


Scheduling 

Lateness 

Shop  fJ.ow  tine 

Number 
of  jobs 

Average 

rule 

Mean    Std. 

Dev. 

Mean    Std.  Dev, 

viumber 

lane 

of  job? 
in  the 
shop 

SH0PN(ERL.4)  17.81  79,99  92.08  53.18 

SHOPN(ERL,S)  1,5  20.34  72.92  48.87 


273     89.01 


195 


72.11 


S/OPN(ERL,S)    365,3     115.2      447,1      127.8 


500    374.63 
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The  interesting  results  are  summarized  below. 

1.  The  graphs  of  inter- arrivals  and  later-departures  show  that 
even  after  600  jobs  have  departed  from  the  ^cio^,    the  curves  are  not 
as  close  as  they  are  in  the  case  of  Experiment  I.   This  is  logical 
since  with  the  increase  of  the  value  of  the  parameter  K  of  Erlang 
distribution,  the  processing  times  are  increased.   Hence  more  jobs 
arrive  in  the  shop  than  those  departed.   Consequently  an  explosion 
could  take  place  in  the  shop. 

2.  The  shortest  irominent  operation  rule  is  more  effective  in 
the  case  of  Erlang  8  than  that  of  Erlang  4. 

3.  In  general,  the  performances  of  the  rules  ace   poorer  than 
those  obtained  when  the  exponential  processing  times  are  assumed. 


aiAPTER  V 
SUMI-lyVHY  MiH   CONCLUSIONS 

The  job-shop  simulation  is  a   process  in  which  different  scheduling 
rules  i-inder  different  experimental  conditions  are  evaluated.   A  simu- 
lation raodel  can  nearly  duplicate  the  conditions  prevailing  in  real 
time  job-shops.   The  purpose  of  this  thesis  is  to  (1)  present  a 
consistent  and  more  comprehensive  review  of  the  most  important  job-shop 
simulation  experin.ents  conducted  since  1960;  (2)  develop  a  method  of 
assigning  job  due-dates,  and  testing  it;  (3)  evaluate  the  efficiency 
of  GASP-II  in  handling  simulation  experiments;  and  (A)  study  the  effect 
of  different  processing  time  distributioa  on  the  performance  of  the 
rules. 

The  design  of  the  simulation  model  is  influenced  by  the  viork  of 
Conway  [12] ,  since  it  is  necessary  to  compare  the  results  on  the  basis 
of  same  experimental  conditions.   However,  limitation  of  computer  tines 
has  been  a  deciding  factor  wliich  influenced  the  experimental  conditions 
in  the  present  investigation.   A  surtimary  of  this  experiment. 3 1  ir.^esti- 
gation,  referred  to  as  Experiment  16,  along  with  the  results  cbtcu'-red 
are  presented  in  the  next  page. 

GASP-II  can  efficiently  handle  job-sliop  probletrs  of  the  size  as 
big  as  15CC  jobs  and  9  machines.   The  average  computer  execution  time 
for  a  simulation  run  of  1300  departed  jobs  (corresponding  to  1500  arrived 
jobs)  is  approxim.ately  20  minutes. 
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EXPERIMENT   16 


INVESTIGATOR  Present   Thesis 

COIIPUTER/PROG.    LANG.  IBM  System   360/50   -   GASF-II 

I'     SHOP    OiARACTERISTICS 

1.  Shop  structure 

2.  Shop   type 

3.  Number  of  groups 

4.  Nmaber  of  machines/group 

5.  Shop  utilization 
II :  EXPERIMENTAL  CONDITIONS 

1.  Distr.  of  job  arrivals 

2.  Distr.  of  job  routings 


Network  of  single  queues 

Job  shop 

9 

1 

=  93% 


3.   D:;'-^t.r.  of  Number  of 
op orations /j ob . 


Exponential  inter-arrivals,  mean  -1,C 

Uniform  [1 ,9] 

Geometric    [4,12]    with  iriean  -   9, 


4.  D.    :r.    of   ^.xocessing  time        Erlanj    (K-1,4,8),   mean  -   0 

5.  Saop  utili nation  1 

6.  Ir.itial  Sliop   State  Empty 

7.  Initial   number    of  - 
released  jobs 

8.  R'j:>in  length  - 

9.  Nurber  of  jobs    considered        Intermediate  statistics   are   col- 


for   collecting   statistics 


10.  R--.r.-out   length 

11.  Jobs/run  cr    simulation 
tiviie 

12.  Due-date  inultjp].ier 


lected  after  every  100  jobs  art; 
departed.   Final  statistics  are 
collected  after  1300  departed  jobs; 
houever^  in  the  case  of  Erlang  4 
and  S  it  is  600  den^rtc-d  jobs. 


1300  and  600  departed  jobs 


1.   Proportional  to  work  content  of  job 
2^   Proportional  to  v/ork  content  and 
number  of  operations/job. 
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EXPERIM:NT  1G  (Contd.) 

III.   Scheduling  rules:   RUL103,  105,  114,  116,  226 
TV.   Measures  of  performance:   NESlOl,  104,  109,  110,  201 
V.   Results: 

1.  The  shortest  imminent  operation  rule  (RUL105)  is  superior 
to  others  in  reducing  job  lateness  and  shop  flow  time. 

2.  The  procedure  in  which  the  due-date  allowance  is 
proportional  to  the  number  of  operations  axid  work  content 
of  the  jobs  has  proved  to  be  beneficial  in  the  case  of  the 
non-due-date  rules  such  as  first-come  first-served  and 
shortest  imminent  operation  tire  rules.   However,  there 

is  no  effect  on  the  performance  of  due-date  rules  (RUL114, 
116,  226). 

3.  The  perfori.ance  of  the  shortest  iraminent  operation  time  and 
slack  per  remaining  num.ber  of  operations  rules  are  affected 
badly  V7ith  the  use  of  the  processing  time  distributions  of 
Erlang  with  parameter  K  equaling  4  and  8.   However,  the 
performances  under  Erlang  8  are  better  th.an  those  obtained 
under  Erlang  4. 

4.  The  GASP-II  v;orks  equally  good  for  large  size  of  shop 
problems.   The  simulation  computer  execution  time  of  runs; 
ranging  from  600  to  1300  departed  jobs,  varies  from  18  to 
23  mi'vutes  depending  on  the  rules  tested. 
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APPENDIX  A 


SCHEDULING  RULES 


Til  is  appendix  liicludes  iiluiost:  a  comprehensive  list  of  scheduling 
rules  along  wiirh  the  corresponding  definitions  and  raathematical 
expressions.   An  effort  has  been  made  to  classify  these  rules  into 
local  static  (LS) ,  local  dynamic  (LD) ,  global  static  (GS) ,  and  global 
dynamic  (GD)  types;  however,  it  is  very  difficult  to  be  conclusive. 
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NOTATION 

j  job    index;      i    -    1,    2,     ...,    J 

j  particular  job 

K,  total  iiuraber  of   operations    for  job  j 

fe  operation   index   for  a  job  j,      fe   =   1,    2,    ...,    K. 

k  particular   operation 

M  total   number   of  macliines 

rn  machine   index;      m   -    1 ,    2,     .,.,    M 

n  sa;Tip].e   si'/,e 

R.(t)  priority  index  of  job  j  at  time  t 

t.i  processing  time  of  job  j  for  operation  k 

D.  due  date  of  job  i 

3 

D_.,  due  date  of  job  j  for  operation  k 

r  ;  ready  time  of  job  i  for  operation  k 


1 


:lme  at  v;hii:h  job  j  arrived  at  shop 


J 
N.,^(t)  nuinber  of  jobs  in  queue  while  job  j  is  ready  for  its  operation 

1  !^ 


k   at    time   t 


0,1  (t)    total  vjork   content   of   a  ([ueue  v/hile  job   j    is   ready   for   its 
operation   'i    it    time   t 

Q.,  (t)    total  expected  vjork   content   of   a   queue  v/li.lle  job   j    is    ready   for 
•.ts    operation  k   at   time   t 
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slack  time  of  job  j 


X.|   random  variable  assigned  to  job  j  for  operation  k 
W.I   exnccts'd  v/alLln'.^  tir.u'  of  i  ob  i  for  oporaticn  fe. 
V.   monetary  value  of  job  j 
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RULIQI:      RANDOM  RULE    (LS) 

Schedule  nexl:    the  job   v/hich   has    the   loi>/est   value   of    a   random   priority 
index.      The   priority   value   are   assigned   randooily   as    each   nev;   job   joins 
the   queue  before  a  machine. 


R.(0    =X.£ 


RUJL102_:   FTRST  IN  SHOP  FIRST  SERVED  (LS) 

Schedule  next  the  job  v;hich  has  the  earliest  arrival  date  at  the  shop. 

The  priority  nurabers  are  given  in  increasing  order  as  the  job  arrive 

in  the  shop.   The  job  vhich  has  the  lowest  priority  number  is  scheduled 

next. 


R.(t)  -  r.. 


RULIOS:   FIRST  IN  QUEUE  FIRST  SERVED  (LS) 

Schedule  next  the  job  which  entered  the  queue  first.   The  priority  numb 
are  given  in  increasing  order  as  each  job  arrives  at  machine.   The  job 
having  lov/est  value  of  priority  number  is  scheduled  first. 


R,(t)  -  r. 


k 
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RUL104:   FIRST  COME  IN  QUF.UE  FIRST  SERVED  WITH  SHORTEST  OPER.\TION 

OVERLOAD  PROTECTION  (LS) 
Schedule  noxt  the  iob  v/hich  ar:;ived  first  in  the  queue  tuilcss  the 
nuTfibc-r  cf.    jobs  waiting  in  the  queue  has  not  crossed  the  value  a. 
If  tha  nvurtbor  of  jobs  in  ihe  queue  is  equal  to  or  greater  than  a, 
the  job  with  the  s'uortest  operation  processing  time  is  scheduled 
first. 


R.(t) 


t.r  ,    otht 


lerwise 


_FaJIA05:   SHORTEST  IM>IINENT  OPERATION  TIME  (LS) 

Schedule  next  the  job  wliich  has  the  shortest  imminent  operation  processing 

ti;;!e.   In  other  words,  the  job  with  short  operation  is  given  first  pref- 


erence. 


Rj(.0=.n.^ 
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RUL106;   IWO-CUJSS  SHORTEST  OPERATION  (LS) 

Schedule  next  the  job  frora  shOi.t  operation  class.   Jobs  being  classified 
into  two  classes:   shoi't  job  class  and  long  job  class.   Preference  is 

given  to  short  job  class.   A  job  is  called  a  sh.-^.rt  job  if  its  inmiinent 

cperatiot-i  processing  time  is  less  than  the  specified  para'.iicter  a  which 

is  mea:  ired  in  units  of  time.   U'ithin  each  class  the  tie  is  resolved 
by  first  come  first  served  rule. 


lyt) 


ill 


b+r  .  r 


^jfe<  ^ 


t .  r   >  a 


RUL107_:   TOREE-CLASS  SHORTEST  OPERATION  (LS) 

Schedule  next  the  job  according  to  the  sliortest  operation  rule  from  a 
set  of  jobs  v;hich  will  arrive  at  the  queue  of  size  equal  to  or  greater 
than  a  constant  a.   All  v/aiting  jobs  are  divided  into  three  classes: 
the  first  class  consists  of  jobs  having  their  next  operation  on  a 
machine  with  a  queue  of  size  equal  to  or  greciter  than  a.   Tlie  second 
class  consists  of  jobs  whose  subsequent  operation  is  to  be  processed 
on  a  machine  with  a  queue  of  size  less  than  b,  where  b  >  a^   The  third 
class  consists  of  the  remaining  jobs.   The  first  clasr^  has  preference 
over  the  second  one  which  has,  in  turn,  preference  over  the  third  class 


R.  (t-) 


'jk 


ik  ' 


^ji+i'^>.-^ 


where  a  and  b  are  positive  integers  such  that  b 


csu 


RUL108;   TRUNCATED  SHORTEST  OPEMTION  (LD) 

Schedule  next  ':lie  job  which  has  the  short'.est  linminent  operation  processing 

time  unless  there  is  a  job  in  the  queue  which  lias  waited  more  than  or 

equal  to  certain  time  units,  a.   U'hen  a  =  0,  this  rule  becomes  same 

as  first  in  queue  first  served  rule. 


ft,,:,    (t-r.p    <    a 

(t)    -  - 
J 


3'^  J 

R. (t)  -  { 


0  ,    (t-r^p  >.  a 


RUL1Q9  •■   LONGEST  IMMINENT  OPERATION  TIME  (LS) 

Schedule  next  the  job  which  has  the  longest  iri-iminent  operation  processing 

time. 

Priority  number  is  given  inversely  proportional  to  the  operation 
processing  time. 


R.(t)  =  1/t.r 


RULllO:   LONGEST  OPERATION  WITH  SHORTEST  OPER,\TIOri  OVERLOAD  PROTECTION  (LS) 
Schedule  next  the  job  with  the  longest  operation  processing  time  unless 
the  number  of  jobs  waiting  in  the  queue  has  not  reached  a  certain  number 
equal  Lo  or  greater  than  a.   If  the  nuni.ber  of  jobs  waiting  is  equal  to 
or  greater  than  a,  schedule  the  job  v/lth  the  shortest  operation  processing 
time. 

f 1/t.r,   N.r(t)  <  a 

R.(t)  -  { 

t  .p      ,        otherv/lse 
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RULlll:   SUBSEQUENT  OPEPJVTION  (CS) 

Sc'p.Gclule  next  the  job  whose  subsequent  operation  is  to  be  processed 

on  a  vacchine  v/i  th  a  critical  queue  ;uid  has  shortest  operation  processing 

time  among  tiie  jobs  which  are  qualifying  for  the  next  critical  queue. 

A  critical  queue  is  the  one  v;ith  less  than  a  certain  time  units  of 

work  V7aiting  to  be  processed,  a.   If  there  is  no  critical  queue, 

select  the  job  in  order  of  their  arrival. 


jk    '  ^J,fe+1 

J 

[r.r  ,    otherwise 


RUL112;   LONGEST  SUBSEQUENT  OPERATION  TIME  (LS) 

Schedule  next  the  job  whose  next  operation  has  the  longest  processing 

time.   Priority  numbers  for  the  current  operations  are  assigned  inversely 

proportional  to  the  subsequent  operation  processing  tirae.   Preference 

is  given  to  jobs  on  their  last  operation  by  assigning  a  very  low 

va.luo  of  priority  index. 


[l/raax[t.r],   k  -  K. 


1'^- 


5/ 


RUL113;   OPERATION  SLACK  TIME  (ES) 

Schedule  next:  the  job  which  has  the  minimum  operation  slack  time. 
Tlie  operation  s.lack  time  is  tiie  difference  between  the  present  and 
previous  operation  due  dates  minus  the  operation  processing  tinie. 
Operation  slack  time  is  the  time  for  the  operation  to  make  up  if 
it  is  already  late. 


R. (t)  -  D.r  -  D.  r  ,  -  t.r 


RUL114:   REMA.IN1NG  JOB  SLACK  TIME  (LD) 

Schedule  next  the  job  v;hich  has  the  minimam  job  allowance  or  free  time 
before  the  job  due  date.  The  job  allowance  is  the  time  left  after  all 
the  operation  processing  times  of  a  job  are  deducted  from  the  due  date 
minus  the  present  time. 


K. 

3 


R   (t)   :=^   D    -   t   "     I.       t   , 


RUL115:   EARLIEST  START  DATE  (L3) 

Schedule  next  the  job  v/hich  has  tlie  earliest  planned  start  date.  A 
planned  start  date  on  v;bich  a  job  should  start  in  order  to  be  completed 
on  time  is  determined  b\'  miultiplying  the  total  number  of  operations 
by  arbitary  allowance. 
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RUL116:   EARI.IEST  JOB  DUE  DATE  (LS) 

Schedule  next  the  job  which  has  the  earliest  job  due  date.   Priority 
numbers  are  assigned  proportional  to  job  due  date.   Jobs  •.•.hich  are 
supposed  to  be  completed  ear]y,  receive  a  highior  priority. 


R.(t)  =  D. 
3  3 


RyUJ_7:      OPERATION   DUE   DATE    (LS) 

Schedule   next    the   job   v.iiich   has    the   earliest    operation   due   date.      Th.e 
operation  due   date   is    computed  as    the  sum  oi;    the  job   arrival   time   at   the 
shop   and   equally  spaced  job   due  date  per  operation. 


R.(0    -r.^+(D.    -r.p    fe/K, 


RUL118;   EARLIEST  OPERATION  DUE  DATE  AND  FROPORTIONM.  SLACK  TIME  (LS) 
Schedule  next  the  job  which  has  the  earliest  operation  due  date.   The 
operation  due  date  is  computed  as  the  sum  of  the  previous  operation 
due  date  and  the  assigned  shop  time  of  job  proportional  to  the  operation 
processing  time. 


R,  (t)  =  D.  r  ,  +  (D,  -  r.J  --— Jl- 


j,fe-l     J     jl   K 


aii 


RUL119:   FEWEST  REMAINING  OPERATIONS  (LS) 

Schedule  next:  the  job  v/hich  has  the  least  number  ol:  operations  still 

to  be  perfornied  before  being  departed  from  the  shop.   The  priority 

nuraber  of  operations  reinaining  includliig  the  present  operation  for 

which  it  is  waiting.   Hence  jobs  on  their  last  operations  are  scheduled 

first. 


R.  (t)  =  K.  -  fe  -f-  1 
J       J 


RUL120:   MOST  REMAINING  OPERATIONS  (LS) 

Schedule  next  the  job  \-/hich  has  the  most  number  of  the  unscheduled  operations 
before  being  departed  from  tlie  shop.   Priority  nuuibers  are  assigned  in- 
versely proportional  to  the  number  of  operations  remaining. 


R,  ( t)  =  1/  (K,  -  k   -f  1) 


RUL121:   GREATEST  TOTAL  WORK  (LS) 

Schedule  next  the  job  which  has  the  most  amount  of  work  on  its  routing, 
Priority  numbers  are  assigned  to  th.e  jobs  on  their  arrival  in  the 
shop.   The  jobs  wh.ich  has  the  maximum  total  operations  processing  time 
is  given  highest  priority.   The  priority  remains  static,  i.e,  does  not 
change  v/ith  time 

K. 

R.(t)  --=  1/y   t.. 
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RUL122;      LEAST  WOilK  REMAINING    (LS) 

Sci\edule   next    the   job  which  has    the   least    .-imouDt    of    remaining  v;ork   to 

be   performed    including    thn   present   operation.      The   amount    of   work 

to  be   perforraed   on   a  job    is    the   suni  of    all   of    its   operations    processiiig 

times. 

K, 

^/t)    -.      l]      t   ^ 
J  k=iz       J'^ 

RUL123 ;   MOST  WORK  REMAINING  (LS) 

Schedule  next  the  job  which  has  the  most  remaining  work  to  be  performed. 
In  other  words  the  job  whose  total  remaining  operation  processing 
times  is  a  maximum,  receives  highest  priority. 

K. 

1 

R.(t)  - 1/  y.  t., 

RUL124;   FEWEST  JOBS  IN  NEXT  QUEUE  (GD) 

Schedule  next  the  job  whose  subsequent  operation  is  going  to  be  processed 

on  a  machine  with  a  queue  having  the  least  number  of  jobs  waiting.   If 

tvv'o  or  more  jobs  qualify,  the  tie  is  resolved  by  the  r.hortest  operation 

ruJ.e. 


P.jCt)  »Nj,fc-,l(0 


ou 


RUL12  5;       LEAST  WORK    EN  NEXT  QUEUE    (GD) 

Schedule   next    the  job  whose  subsequent   operation   is   going    to  be 
processed   on  a  machine  with   a  queue  having    the  ininjmum  sum  of 
processing    times   of   iill  vvaiting   jobs. 


V'^-'^iM^'^ 


RUL126:   LEAST  EXPECTED  WORK  IN  NEXT  QUEUE  (GD) 

Schedule  next  the  job  whose  subsequent  operation  is  going  to  be  per- 
forraed  on  the  machine  of  which  the  queue  hns  the  minimum  sum.  of  v7ork 
waiting  and  work  expected  to  arrive  soon.   Expected  work  ia  the  sum 
of  operations  processing  time  of  all  those  jobs  which  at  present  being 
processed  on  different  machines  and  will  arrive  soon  in  the  queue 
under  consideration. 


V^  ^^j,^.-i^^^  -^^j,i;+i^^> 


RUL127;   FIRST  HALF  PREFERENCE  (LS) 

Schedule  next  the  job  which  is  yet  to  be  processed  on  more  than  half 
its  total  number  of  operations.   This  rule  divides  the  jobs  which  arc 
waiting  for  processing  into  '  /o  classes  depending  ou  the  correspox;dlng 
number  of  operations.   A  tie,  if  exists,  is  broken  by  first  in  queue 
first  served  rule. 


fl   ,    (K.  -  'm-D  <  K./2 
J        -  .1 
R.(i:)  =-  -i 
J       ! 

(0   ,    (K.  -  f'.-^l)  >  K./2 
J  J 
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RUL128;      LAST  HALF  PKEFERENCE    (LS) 

Schedule  next    the  job  which   is   yet    to  be   processed  on   less    than,  or 

equal    to  half   of    its    total  nuinber   of   operatior^s .      If    a   tie   exists, 

schedule    tVie   j  cb   t-.'hich   iirrlved    ear""  ier   in    f'he   qn-Gae,      If    there   is 

nc   job    to   qual:ify    ti\c   above    conditioa,    sclicdiile    the   jobs    i-iccordir.g    to 

first    in   queue   first   served    rule. 


'O,         (K.    -   k  +   1)    <_  K../2 
R.  (t)    - 

[1,         (K.    -   fe  ^-l)      >    K./2 


faJ_IJL_2_9:   LATE  OPEIUTION  PREFERENCE  WITH  SHORTEST  IMMINENT  OPEP^\TION 

TIME  (LS) 
Schedule  next  the  job  which  is  behind  operation  due  date.   If  there  is 
no  such  job  in  the  queue  then  the  job  which  has  the  shortest  immineat 
operation  processing  time  is  selected  next. 


t.r  ,    r.,:  <  r.^  -h  (0.  -  r.)  k/K. 


R.  (t^ 
I 


[0    ;,    other^^ise. 
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RUL.130;   LAST  OPEMTION  OF  LATE  JOBS  PREFERENCE  WITH  SHORTEST  J.I4MINENT 

OPERATION  TIME  (LS) 
Schedule  next  the  job  which  has  the  stiortesl:  linminent  operation  processing 
time  unless  there  is  a  late  job  waiting  in  the  queue  for  its  last 
operation  ^-zhose  lateness  can  be  avoided  by  scheduling  it  first. 


R.(t)  - 


U.   ,    G  <  P   for  fe  -  K.  and  U.  >  0 


G  >  P 


where. 


u.  -  D.  --  t  -  t.r 


P   =  nun  [t^] 
J 


G  =  min  [U.]  for  all  jobs  with  fc  ^=  K.  and  U.  >  0 
J  3      J 


RUL201:   LAST  OPERATION  PREFERENCE  WITH  POSITIVE  OPERATION  SLACK  TIME  (LS) 
Schedule  next  the  job  which  has  the  minimum  sum  of  weighted  remaining  job 
slack  time  ard  weiglited  imminent  operation  processing  time.   The  weighting 
factor  for  processing  time  is  the  function  of  work  content  of  the  queue. 
The  weighting  factor  of  job  slack  time  gives  preference  to  the  last 
operation  of  the  job. 

R.(t)  -  a'^  (D.r  -  D.  r  T  -  t.r)  +  (b  y  t.,)  t.r 

where 

a  >  1 


(  0,  k   i   K. 

^-l,    \z   -  K.   and   (D .  P  -  D.  7    ,  -  t.T)  >  0 

I  1,    fe  -  K.   and   (D.P  -■  B .  7    ,  -  t.D  <  0 
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RUL202:      LAST  OPERATION  PREFERENCE  WIT^I   JOB   SLACK  TIME   FACTOR    (LS) 
Schedule   next    the  job   with    the   siriallerjt    suv.\   of    inimin.^nt   operation 
slack    time   and  weighted  work   in    tl^e   queue  rainus    a  value,    v   which 
takes   a  positive  nuraber   iC    the   current   operation   is    the   last   operation 
and    its    operation  slack   time   is    less    than  u,    other. .'iye  v  has   value 
equal   to   1. 


«j^')  =  (°jfe-Oj,fe-i-^jP  +  (^I'jP  'jfe-^ 


where 


V   >    0    ,         k   =--  K.    and    (D.P    -   D.    T    ,    -    t.f)    <   u 


V   =   0    ,        f.i   ^  K.    or      (D.r   -  D.    P        -   t.p    >   u 
J  j'i  1,'^--!-  3'^     - 


RUL203 ;   LAST  OPERATION  PPvEFERENCE  WITH  OPERtVTICN  SLACK  TIIiE  AND  PROCESSING 

TIME  FACTORS  (LS) 
Schedule  next  the  job  which  has  the  minimum  sum  of  imminent  operation 
slack  time  and  weighted  imminent  operation  processing  time.   The 
weighting  factor  is  a  function  of  work  content  of  the  queue  and  assigns 
a  higher  priority  to  a  job  v/hich  is  waiting  for  its  last  operation. 


L(t)  -  (D.r  --  D.  r  .  -  t.p  -I-  [xb  y  t.r]  t.r 

J  J'^  l,f^-l  3'^  7,      jfe        J'^ 


wner 


r 

!0.        k   ■-  K. 

i'  "1 

1  1  •  ^ 

1 ,        ocherv.'ise 
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RUL204:   LAST  OPERATION  PllEFEPvENCE  WITH  DIFFERENCE  BET^^?EEN  PROCESSING 

TIME  AND  SUBSEQUENT  PROCESSING  TIME  (LS) 
Schedule  next  the  job  which  is  waiting  for  its  ].ast  operation.   If  no 
such  job  is  i/aitiiig,  select  the  job  with  the  least  difference  between 
its  weighted  current  processing  time  iind  v/cighted  subsequeat  operation 
processing  time  wliere  a  is  a  weighting  factor.   When  a  =  1,  this  rule 
becomes  the  shortest  operation  rule;  however,  when  a  =  0,  it  becomes 
the  longest  subsequent  operation  rule. 

lat.r  -  (1-a)  max  [t.,],   fe  =  K. 

RUL203;   SHORTEST  OPERATION  WITH  PREFERENCE  TO  ARRIVING  JOSS  WITH 

SHORTER  PROCESSING  TIME  (GS) 
Schedule  next  the  job  with  the  shortest  operation  processing  time  with 
a  'hold'  if  there  are  few  jobs  in  the  queue  and  a  job  with  shorter 
processing  time  will  arrive  soon. 


RAt)  -  t. 


J 


k 


No   assigameiit   is  made,    i.e.    the  machine   is   held   idle,    if   Q.T    <    a  and   if 

J'-  " 

there  is  a  job  j  nov7  being  processed  for  its  (fe  -  1)  operation  that  will 
come  to  this  queue  for  its  K-th  operation  with: 


^ji-1  '   "^!^  f^jfe^  -  ^ 


where  b  is  a  constant. 
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RUL20_6:      DIFFERENCE  BETWEEN   AVAILABLE   TIME  UNTIL  DUE    DATE   MB   I>IMINENT 

OPERATION   TIolK    (LS) 
Sched'jle   next    the   job  with    l:b.e  niiniiaum   dif Id'eiice   of   weighted    reiiiaining 
flow   tiiTie    (job   due   date   minas    present    tine)    aiid    imnilreni:    opi^ration   time, 

R.(t)    -  w(D.   -   t)    -  t.r 

where  w  is  waiting  factor. 

RUL207_:   ALTERNATE  OPERATION  HEURISTIC  (LS) 

Schedule  next  the  job  which  becomes  critical  as  a  result  of  scheduling 
another  job  according  to  a  specific  priority  rule.   If  processing  of 
til'.-:  critical  job  first  create  another  critical  job  in  the  queue, 
schedule  next  a  j ob  according  to  the  specific  priority  rule.   A 
critical  job  is  the  one  which  is  already  late  or  soon  will  be  late. 

RUL208:   LOOK  AHEAD  HEURISTIC  (GS) 

Schedule  next  the  job  v/hich  is  critical  (late  or  near  late)  and  is  due 
to  join  the  present  queue  at  some  future  time  before  the  completion  of 
the  oneration  of  the  job  which  has  been,  selected  by  a  priority  rule. 
If  the  processing  of  this  future  critical  job  crccite  another  critical 
jobs  in  the  present  queue,  apply  the  current  priority  rule  and  neglect 
f u t u r e  cr  1 1 1  c a  1  job. 


9?. 


RUL209 ;   INSERT  HEURISTIC  (LS) 

Schedule  next  the  job  v^hose  operati.'n  processia.'.3  time  can  Tit  into 
the  idle  time  created  by  look  ahead  rule.   Due  to  look  ahead  rule 
the  jobs  at  the  present  machine  are  withheld  from  being  processed, 
since  critical  job  is  due  to  arrive  at  the  queue  soon.   Tliis  causes 
an  idle  time  of  the  corresponding  machine. 

RUL210:   TIME-TRANSCENDING  SCHEDULE  HEURISTIC  (LD) 

Schedule  next  the  job  v/ith  the  top  priority.   Re-evaluate  the  priorities 

every  time  a  job  joins  the  queue. 

RUL211:   SUBSET  OF  CRITICAL  JOBS  HEURISTIC  (LS) 

Schedule  next  the  job  from  a  subset  of  critical  jobs  according  to  a 
priority  rule  or  set  of  rules,  then  schedule  the  raiiaining  jobs 
around  them.   In  this  rule  schedule  is  laid  out  by  jobs  instead  of 
job  operations. 

RUL212:   RE-DO  WITH  ADJUSTED  DUE  DATES  HEURISTIC  (LD) 

Schedule  next  the  job  according  to  the  priority  rule.   When  the  schedule 
is  completed  check  for  late  jobs,  decrease  their  due  dates  by  the  time 
the  job  v/as  late  and  schedule  again.   This  will  give  chance  to  the 
previously  J.iite  jobs  to  be  completed  in  time. 

_RUL213:   M^ilN'IPULATION  HEURISTIC  (LD) 

Schedule  next  the  job  as  determined  by  the  Oantt  chart  until  a  better 
schedule  is  obtained.   'Hiis  heuristic  is  similar  in  nature  to  the 
heuristic  of  re-do  with  adjusted  due  dates. 
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RUL214:      ONE   CLASS   MONETARY   VA].VE   OF   TLIE  JOB    (LS) 

Schedule  next    the  job  with   the  highest  inonatary  value  among   jobs  waiting 

in    the   queue.      Priority   index   is   inversely   proportional   to    the  valvie 


;f   job. 


I^.  (t)    -    l/^. 
J  J 


RUL215;      TWO- CLASS   MONETARY   VALUE   OF   Tlffi   JOB    (LS) 
Schedule   n?xt   the  job   from  higher  monetary   value   cl^iss.      Jobs    in 
queue   are    divided   into    tv/o    classes   depending    on    their  raonetary   value. 
All  jobs  with  higher   raonetary  value  have  higher   priority.      Within 
each   class,    the  job  which   arrived   first   in  queue   is   selected   first. 

! t . r   +  max    [  t ,  r ]  ,        V .    <    a 

where  a  is  a  specified  parameter  measured  in  monetary  value. 

RUL216;   TtTREE-CLASS  MONETARY  VALUE  OF  THE  JOB  (LS) 

Schedule  next  the  job  which  arrived  first  in  the  highest  monetary  value 
class.   Jobs  are  classified  into  three  classes:   high,  medium  and  low 
depending  on  their  raonetary  values.   V/ithin  each  class  jobs  are  selected 
v;.i  fch  first  come  first  served  rule.   In  other  words  jobs  arc  selected 
from  th.e  highest  monetary/  value  queue  first. 

R.  (;:)  -  "(pl-r.r,      b  <  V.  <    a 
J      \        '-^'^  ~  J  ■  ' 

•^q-l-r.r      V  <  b 
J '-       J 

where  a  and  b  are  specified  parameters  measured  in  monetary  values,  a  >  b, 
aiid  p  and  q  are  constants  such  uliaL  q  '->  p. 
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R_UL2_17_:   E:VRLIESX  DUE  DATE  AND  SUBSEQUENT  QUEUE  SIZE  (CD) 

Scliedule  next  the  job  which  has  the  earliest  job  due  date  and  iv'hose 

subsequent  operation  is  going  to  be  processed  on   a  raachiac  with  a 

queue  having  less  than  a  certain  number  of  v/aiting  jobs,  a.   If  such 

a  job  does  not  exists,  apply  the  simple  earliest  due  date  rule.   This 

compound  rule  divides  waitiiag  jobs  into  two  groups:   one  consisting  of 

jobs  whose  subsequent  operation  is  to  be  performed  on  a  machine  with 

queue  of  a  size  lass  than  a.   The  second  group  consists  of  jobs  v.^bose 

subsequent  operation  is  to  be  processed  on  a  machine  with  a  queue 

of  a  size  greater  than  a.   The  jobs  in  the  first  group  is  given  preference 

over  those  in  the  second  group.   Furthermore,  the  jobs  in  Loth  groups 

are  scheduled  according  to  the  earliest  due  date  rule. 


R.(t)  - 
J 


D. ,     fc  -  K.  or  N.  r , , (t)  <  a 


D,  +   n,  otherwise 
(  1 


RUL218:   FEWEST  JOBS  IN  NEXT  QUEUE  AND  SHORTEST  OPEMTION  (GD) 
Schedule  next  the  job  v/hose  subsequent  operation  is  going  to  be 
processed  on  a  miachine  with  a  queue  having  a  number  of  jobs  less  than 
a,  i.ti.,  a  short  queue  is  available.   If  all  the  queues  have  a  number 
of  jobs  more  than  a,  i.e.,  no  short  queue  is  available,  then  select 
tiiC  job  by  applying  the  shortest  operation  rule. 


R.(t)    -^     \ 
J  I 


J    7  -^   '    -^ 


■■\t. 


""^Mi^'^  "  ' 


otherxv'ise 
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RUL219:   REMAINING  JOB  SLACK  TIME  MD  SHORTEST  OPEMTIOM  WITH  CONTROL 

OVER  JOB  DELAYS  IN  SHOP  (LS) 
occKidule  novL  Lire  job  with  the.  shortest  operation,  processing  time  from 
the  set  of  jobs  having  F  <  0  where: 

K. 
F=  {D   -  O  -     f    t  ,  -  c 

■^  fe  =  fl    ^ 

and  c  is  a  control  parameter.   If  no  such  job  is  v/aiting,  simply  schedule 
a  job  according  to  the  shortest  operation  rule. 


R.Ct) 
J 


^jfe  ■      '^° 


I t . r  ,     otherwise 
^  3  re 


RlJL_22_q:   MODIFIED  REMAINING  JOB  SLACK  TIME  AND  SHORTEST  OPERATION  WITH 

CONTROL  OVER  DELAYS  (LS) 
Schedule  next  the  job  with  the  minimum  value  of  F  from  set  of  waiting 
jobs  having  F  <_  0  where: 

K. 

F  -  (D   -  t)  -  I.    t  ,  -  c 

and  c  is  a  control  parameter.   If  xio  such  a  job  is  waiting,  schedule 
a  job  according  to  the  shortest  operation  rule. 

F,      F  <  0 


R.(t)  H 


t  .  r  ,  otherx/ise 

-5  ,' ' 


96 


RUL221;   ALTERNATING  SHORTEST  OPERATION  AND  FIRST  COME  FIRST  SERVED  (I.D) 
Schedule  next  the  job  according  to  shortesc  aparation  rule  for  certain 
fixed  period  of  time,  a,  and  then  schedule  those  accoiding  to  first 
come  first  st-fved  rule. 


[r.r  ,      for  certain  period  of  time,  a 


R  (t)  H 


'jfe 


othert.'/ise 


RUL22_2:   RATIO  OF  PROCESSING  TIME  TO  NL^IBER  OF  OPERATIONS  REMAINING  (LS) 
Schedule  next  the  job  which  has  the  miniraum  ratio  of  its  processing 
tiiiie  to  the  weighted  number  of  operatioi\s  remaining  including  the  current 
operation.   Hence  for  the  same  operation  processing  tirae,  the  job 
which  has  Kore  regaining  operations  is  given  higher  priority.   Ivhea  a, 
the  weighting  factor  is  such  that  a  ->  «^ ,  this  rule  becomes  the  most 
operation  rule;  however,  when  a  -  0.  this  rule  becomes  shortest  operation 
rule. 


R,(t)  -  t.r  /(K.  -  k-M) 


a 


RUL223;   RATIO  OF  PROCESSING  TIME  TO  WORK  REMAINING  (LS) 
Schedule  next  the  job  which  has  the  minimum  ratio  of  processing  time  to 
work  remaining  raised  to  power  a,  where  a  is  a  v;eighting  factor.  \Tviea 
a  =  0,  chis  rule  becomes  shortest  operation  rule;  however,  v/heri  a  -^  °= , 
this  rule  becomes  the  most  v/ork  remaining  rule. 

K. 


^(O  -  t  r  /(  I.   t^) 


k=k 
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RUL224;   RATIO  0?  PROCESSING  TIME  TO  TOTAL  JOB  WORK  (LS) 
Schedule  next  the  job  v/hich  has  Lhe  snallest  ratio  of  processing  Lime 
to  the  total  job  v/oc-.   Total  work  of  a  job  ±s    the  sum  of  all 
operations  processing  time. 

K. 

R.(t)  -  t.,^  /      )       t., 
J       2k        ,,^1  Jfe 

RUL22.5:   KATIO  OF  S'JM  OF  PROCESSING  TIME  AND  WOM-C  IN  NEXT  QUEUE  TO 

PROCESSING  TIME  OF  NEXT  OPERATION  (GD) 
Schedule  next  the  job  which  has  the  miairaura  ratio  of  sura  of  its  weighted 
processing  time  and  v/eighted  work  content  of  the  next  queue  to 
processing  time  of  its  next  operation.   The  v/ork  content  of  the  next 
queue  is  the  sum  of  the  processing  times  of  jobs  v/aiting  in  queue. 


-      (t   -  ,,) 


RUL226:   RATIO  OF  REMAINING  JOB  SLACK  TIME  TO  REMAINING  NUMBER  OF 

OPERATION  (LS) 
Schedule  next  the  job  which  has  the  smallest  ratio  of  remaining  job 
slack  ?:■  ■■e  ;  (^  the  remaining  operatiotis  including  the  one  for  which 
the  job  is  waiting.   The  remaining  slack  time  is  computed  as  the 
difference  of  job  due  date  and  the  sum  of  remaining  operation  processiag 
time  lainus  the  prese;^^  f  ivno. 

K. 


R.,(t)  -  (D.  -  t  -  h    )l    {^,    -   k+1) 


9  8 


RUL227_:   RATIO  OF  TIME  ELAPSED  SINCE  COMPLETION  OF  PREVIOUS  OPERATION 

TO  REMAINING  NUMBER  OF  OPERATIONS  (LS) 
Schedule  next  the  job  which  has  maximum  ratio  oL"  the  tinae  elapsed  since 
coraoletion  oi  previous  operation  minus  a  v.'eighting  factor  a,  to 
remaining  number  of  operations.   The  v/eighting  factor  a  is  a  function 
of  job  flow  allowance. 

(t  -  r.  p  -  a 

R.(t)  ^- ^'- -■ 

J        (K.  -  fe  +  I)'' 

where  a  and  b  are  positive  integers 

RUL228;   RATIO  OF  REMAINING  JOB  SLACK  TIME  TO  REMAINING  TIME  UNTIL 

JOB  DUE  DATE  (LS) 

Schedule  next  the  job  which  has  the  smallest  ratio  of  free  time  available 

before  the  job  due  date  and  time  remaining  until  job  due  date. 

K. 
R,(t)  .  (D.  -  .  -  J^.p  /   (D.  -t) 
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RUL_22^9:   RATIO  OF  MODIFIED  JOB  SLACK  TIME  TO  REMAINING  TIME  UNTIL  JOB 

DUE  DATE  (LS) 
Sch'-'diile  n.ixt  the  job  which  has  the  minimum  rat  i.o  of  free  tirae  available 
before  due  date,  after  expected  delay  time  for  the  rc;nainiag  operations 
have  been  deducted,  to  the  time  remaining  until  job  due  date. 

K. 
j^-^  (D.  -t) 


where 


t!/  delay  tirae  of  i ob  i  for  operation  k 


RUL230;   RATIO  OF  COST  TO  PROCESSING  TIME  (COVERT  RULE)  (LS) 
Schedule  next  the  job  with  the  largest  ratio  of  delay  cost  to  Imniineat 
operation  processing  tirae.   It  is  assumed  thiit  the  cost  of  delay,  c 
is  simply  the  incrc^iental  change  in  tardiness  of  the  job  order. 


where 


K. 

J 

y  w.(  -  s. 


K. 
3 
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RUL231:   SUM  OF  PROCESSING  TIME  AND  REMAINING  Sf.ACK  TIME  PER  RliMAINING 

OPERATIONS  (LS) 
Sclieclule  next  the  job  which  has  the  minimum  sum  of  weighted  processing 
time  and  weighted  remaining  slack  time  per  remaining  operations. 
Weighting  factor  a,  is  a  function  of  the  number  of  operations  remaining 

K. 

.J 
D.  -  t  --   ^.  t. 

R  (t)  -^   ct  r  +  (1-c) :r-^^\- 

J        J^         (K.  -  fe  +  i)^ 

where 

c  =  a  +  min  [(l-a),  (K.  -  k   H-  l)/b] 

a  and  b  are  positive  integers. 

RUL231=   SUM  OF  PROCESSING  TIME  AND  REMAINING  SLACK  TIME  PER  REMAINING 

OPERATION  (LS) 
Schedule  next  the  job  having  minim.um  ratio  of  weighted  processing  time 
and  weighted  remaining  slack  time  per  remaining  operations  including 
the  one  for  I'^hich  job  is  waiting. 


D. 
J 
R,(t)  -  a  t.r  +  (1-a)  — 
■^  ^'^  (K.  -  k  +   1) 


v/here   a  is    a  weighting   factor. 
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RUL23_3:   SUM  01'  PROCESSING  TIME  MID   WORK  RF,MAINJNG  (LS) 
Schedule  next  the  job  which  has  the  rnininura  surn  of  weighted  processing 
time  and  weighted  work  remaining.   When  a  ==  1,  where  a  is  the  weighting 
f Victor,  this  rule  becomes  the  shortest  operation  processing  time  rule; 
hov;ever,  when  a  -  0,  this  rule,  becoraes  the  least  work  remaining  rule. 


K. 


RUL234:   SUM  OF  PROCESSING  TIME  AND  NWBER  OF  JOBS  IN  THE  SUBSEQUENT 

QUEUE  (CD) 
Schedule  next  the  job  which  has  the  minimum  sum  of  xveighted  processing 
time  and  weighted  number  of  jobs  in  the  next  queue.   When  a  ==  1 ,  where 
a  is  a  weighting  factor,  this  rule  becomes  the  shortest  imminent 
operation  processing  time  rule,  and  when  a  =  0,  the  job  which  will  bo 
processed  for  its  subsequent  operation  at  the  smallest  queue  size  Ls 
scheduled  next. 


R,(0  -=a..,;+  (1-a)  N._£^,(t) 
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RUL235:   SUM  Of  PROCESSING  Tim    ''V.r)   WORK  TM  Ni"XT  QUEUE  (GD) 
Schedule  next  the  job  which  hns  the  snialleiiU  sum  of  weighted  pi-ocessing 
time  and  weighted  work  in  the  next  queue.   Work  In  the  next  queue  i.s 
the  sum  of  operation  processing  time  of  jobs  v;aiting  in  thiit  queue. 


R.Ct)  =  at.,:  -!-  (i-^)  Q.  Mi('') 


RUL236:   SU>1  OF  PROCESSING  TIME  AND  WEIGHTED  EXPECTED  WORK  IN 

SUBSEQUENT  QUEUE  (GD) 
Schedule  next  the  job  \/hich  has  the  smallest  sun  of  its  weighted  processing 
time  and  weighted  expected  work  in  the  subsequent  queue.   The  expected 
work  is  the  v;ork  content  of  the  jobs  \7aiting  and  those  v;hich  will  arrive 
soon  in  the  queue.   The  weighting  factors  being  coinplement  of  each 
other. 


R  (t)  -  at  r  +  (l-a)[Q   T   (t)  +  Q  Z.A':)] 

3  j£  J  ,fc+l  J  ,fc+l 


RUL237;   SUl-I  OF  OPERATION  SLACK  TIME,  PROCESSING  TIME  AND  WOP^  CONTENT 

IN  THE  SUBSEQUENT  QUEUE  WITH  PREFERENCE  FOR  LAST  OPERATION  (GD) 
Schedule  next  the  job  with  the  smallest  sum  of  imminent  operation  slack 
time,  v/eighted  inniiinent  operation  time  and  weighted  work  content  in  the 
subsequent  queue  minus  the  last  operation  indicator,  y  which  takes  value 
greater  than  zero  if  the  job  is  waiting  for  its  last  operation. 

where, 

y  -  0,    k   ^-  K^ 


3 

0<y<c    f^  =  K. 
3 
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RUI^233:   SUM  01'  PROCESSING  TIKE  AILO  WORK  RE,MAINING  (LS) 
Schedule  next  the  job  which  has  the  minimum  sum  of  weighted  processing 
tiwe  and  weighted  work  remaining.   When  a  =  1,  v7here  a  is  the  weighting 
factor,  this  rule  becomes  the  shortest  operation  processing  time  rule; 
hov:ever ,  when  a  =  0,  this  rule  becomes  the  least  work  remaining  rule. 


K. 
R.(r.)=at.^.(lS.J  (i-a) 


RUJ^34:   SUil  OF  PROCESSING  TIME  AND  NUI^BER  OF  JOBS  IN  THE  SUBSEQUENT 

QUEUE  (GD) 
Schedule  next  the  job  which  has  the  minimum  sum  of  weighted  processing 
time  and  weighted  number  of  jobs  in  the  next  queue.   When  a  =  1,  where 
a  is  a  weighting  factor,  this  rule  becomes  Che  shortest  imminent 
operation  processing  time  rule,  and  v/hen  a  -  0,  the  job  v.diich  will  be 
processed  for  its  subsequent  operation  at  the  smallest  queue  size  is 
scheduled  next. 


R.(0  =at.£+  (l-a)  N.__^^j(t) 
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RUL.235:   SUM  OF  PROCESSING  TI?1>:  -'V-T)   WORK  T^i  ^^''XT  QUEUE  (GD) 
Schedule  next    the  job  v;h,lch  hr.s  the  sniallesu  sum  of  weighted  processing 
time  and  weighted  work  in  the  next  queue.   Work  In  the  next  queue  is 
the  sum  of  operation  processing  time  of  jobs  v;aitLng  in  that  queue. 


R.ft)  =  at.,;  -I-  (i-a)  Q.  u.AO 


RUL236:   SUM  OF  PROCESSING  TIME  AND  WEIGHTED  EXPECTED  WORK  IN 

SUBSEQUENT  QUEUE  (GD) 
Schedule  next  the  job  which  has  the  smallest  sum  of  its  weighted  processing 
time  and  weighted  expected  work  in  the  subsequent  queue.   The  expected 
work  Is  the  v7ork  content  of  the  jobs  x/aiting  and  those  which  will  arrive 
soon  in  the  queue.   The  weighting  factors  being  complement  of  each 
other. 


R 


.(t)  =  at.^+  (i-a)[Qj_,:,,,(t)  ■'•Qj_,:^.i(oi 


RUL237;   SUM  OF  OFERiVTION  SLACK  TIME,  PROCESSING  TIME  AND  WOP^K  CONTENT 

IN  THE  SUBSEQUENT  QUEUE  WITH  PREFERENCE  FOR  LAST  OPERATION  (GD) 
Schedule  next  the  job  with  the  smallest  sum  of  imm.inent  operation  slack 
time,  v/eighted  imminent  operation  time  and  v;eighted  work  content  in  the 
subsequent  queue  minus  the  last  operation  indicator,  y  xvhich  takes  value 
greater  than  zero  if  the  job  is  waiting  for  its  last  operation. 

v/here, 

y  -  0,   k  i   Kj 

0<y<c    \l   =   K. 
J 
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RUL238:   SUM  OF  IILMINhNT  OPEIUTION  SLACK  T.1ME  MD  WEIGHTED  OPERiVTION 

PROCESSING  TIME  (LS) 
Scheoule  next,  the  job  which  has  the  smallest  sun  of  in'mincnt  operi\tlon 
slack  tine  anJ  weighted  imminent  operation  processing  time.   The 
weighting  factor  is  the  function  of  the  Vvrork  content  in  the  present 
queue . 

where  b  is  a  positive  integer 

RUL239 :   SLTM  OF  OPERATION  SLACK  TI>IE ,  PROCESSING  TINE  AND  WOPJK  CONTENT 

IN  THE  SUBSEQUENT  QUEUE  (GD) 
Schedule  next  the  job  which  has  the  mj.nitnuni  sum  of  irarninent  operation 
slack  time,  weighted  irarninent,  operation  processing  time  and  v;eighted 
\/ork  consent  in  the  subsequent  queue.   The  weighting  factor  of  operation 
processing  time  is  the  function  of  the  work  content  of  tlie  present 
queue 

R.(t)  -  (D.r  -  D.  r  ,  -  t.r)  ^-  (b  y  t.r)  t.;;  i  cq.  r,.(t) 

where  b  and  c  are  positive  integers. 


104 


RUL240:   WEIGHTED  OBJECTIVE  RULE 

Schedule  next  the.  job  with  maxiruMn  value  ci  nr.  iority  ijitlcx  as  calculated 

by 

R  R 

R.(t)  -   )'   a.  R*.,     y   a.  -  1  and  R- .  <  1 

J      ':i  J  i-j    Jl:      J 


where 


r-1  ^   ^-J    r-1   -^  '^J 


R  -   number  of  rules  combined 
a.  =  weighting  factor  of  job  j 

R" .  ~-  ratio  of  the  'best'  value  for  the  rule  r  to  the  value 
1^3 


of  job  j  for  rule  r. 


Hence,  when  the  rule  considered  is,  say,  shortest  process  lag  ti^ne,  the  job 
with  minir.'.uin  processing  time  v/ill  have  ratio  of  R* .    -   1.00  and  all  other 
jobs  will  have  R* .    <_  1.0  deper.ding  on  the  relation  of  their  processing  time 
to  the  processing  time  of  the  'best'  job. 
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APPENDIX  B 

This  appendix  Includes  the  graphs  of  iutcr-arrivals  and  intcr- 
depar Lures  for  ail  the  experiments  conducted  in  this  investigation. 
Figures  B-1  through  3-5,  B-6  through  B-10,  and  B-.11  through  B-14  correspond 
to  CxperiiTicnts  I,  II,  and  III.  respectively. 

In  addition  to  the  above,  a  computer  output  of  sample  rosult.s  is 
included,   Tliese  results  belong  to  the  shortest  imminent  operation 
rule  tested  in  Experiment  JT. . 
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**GASP  SUMMARY  REPORT** 


**GENcKArED  DATA** 
CODE      GBS.    MEAN     STD.  DEV,    MIN.     MAX, 

STATISTICS  FOR  JPB  LATENESS 


1 

1300 

-35.48 

24.08 

-106.24 

-4.17 

2 

100 

-35.04 

22.97 

-82.80 

-4.17 

3 

100 

-34.6G 

24.48 

-106.24 

-5.11 

4 

100 

-34.10 

23.77 

-78.85 

-4.57 

5 

ICO 

-36.39 

25.22 

-86.09 

-4.79 

6 

100 

-37.60 

25.81 

-93. 7C 

-5.75 

7 

100 

-35. 7C 

24.24 

-85.07 

-5.69 

8 

100 

-37.58 

25.76 

-85.91 

-5.49 

9 

ICO 

-32.52 

21.84 

-73.96 

-5.37 

10 

100 

-33.73 

23.76 

-80.57 

-5.63 

11 

100 

-38.22 

22.39 

-82.76 

-5.40 

12 

ICO 

-37.16 

26.21 

-98.64 

-5.19 

13 

IOC 

-34.09 

23.21 

-89. IC 

-4.88 

14 

100 

-34.53 

2  3.60 

-8  7.03 

-5.67 

15 

MC  VALUES 

RECORDED 

STATISTICS  FOR  SHOP  FLOW  TIME 


16 

1300 

7.54 

5.49 

0.79 

27.41 

17 

ICO 

7.15 

5.14 

0.79 

19.88 

18 

100 

7.16 

5.38 

0.83 

21.08 

19 

100 

7.63 

5.43 

1.06 

23.36 

20 

ICO 

7.66 

5.61 

0.79 

2/. 28 

21 

ICO 

8.56 

6.14 

0.79 

23.12 

22 

100 

7.85 

5.84 

0.79 

22.29 

23 

100 

7.48 

5.23 

0.96 

19.93 

24 

100 

7.48 

5.90 

0.79 

27.41 

25 

100 

6.35 

4.68 

0.87 

22.99 

26 

100 

8.84 

5.58 

0.79 

25.92 

27 

ICO 

8.29 

5.79 

0.79 

20.21 

28 

100 

6.94 

5.36 

0.80 

25.10 

29 

100 

6.61 

4.81 

0.79 

19.18 

30 

NO  VALUES 

RECORDED 

31 

NO  VALUES 

RECORDED 

■*GASP    SUMMARY    REPORT** 
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CODE 


♦♦GENERATED  DATA** 
DBS.    MEAN     STD,  DEV,    MIN, 


MAX, 


STATISTICS  FOR  INTER-DEPARTURE 


32 

13C0 

0.89 

0.81 

0.0 

5.51 

33 

100 

1.04 

0.93 

CO 

4,06 

34 

100 

0.88 

0.77 

0.01 

3,46 

35 

100 

0.81 

0.7D 

0.0 

4.19 

36 

ICO 

0.91 

0.84 

0.0 

4.39 

37 

ICO 

0.84 

0.73 

0.0 

2.82 

38 

100 

0.86 

0.91 

0.0 

5.51 

39 

ICO 

0.87 

0.78 

0.0 

3.81 

AO 

100 

0.79 

0.69 

0.02 

3.44 

41 

NO  VALUES 

RECORDED 

STATISTICS  FOR  INTER-ARRIVAL 


42 

1306 

0.89 

0.42 

0.50 

1.50 

43 

110 

0.95 

0.43 

0.50 

1.50 

44 

1C3 

0.85 

0.41 

0.50 

1.50 

4  5 

96 

0.85 

0.40 

0.50 

1.50 

46 

101 

0.90 

0.42 

0.50 

1.50 

47 

99 

0.86 

C.39 

0.50 

1.50 

48 

102 

0.85 

0.42 

0.50 

1.50 

49 

101 

0.86 

0.44 

0.50 

1.50 

50 

97 

0.81 

0.38 

0.50 

1.50 

**TIVE  GENHRATbD  DATA** 
CnCE    MEAN    STU.  DFV.     TLOB    MAX.    TOTAL  TIME 
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1 

C.92 

0.27 

1155.45 

1.00 

1155.450 

2 

0.94 

0.23 

1156.99 

I.Cj 

1156.990 

3 

0.96 

0.2C 

1154.67 

l.Ou 

1154.670 

4 

0.95 

0.22 

1157.00 

1.00 

1157.000 

5 

0.93 

0.26 

1156.14 

l.CO 

1156.140 

6 

0.95 

0.21 

1156.34 

1.00 

1156.340 

7 

0.94 

0.24 

1146.56 

1.00 

1146.560 

8 

0.92 

0.28 

1153.66 

1.00 

1153.660 

9 

0.92 

0.27 

1156.85 

1.00 

1156,850 

10 

9.47 

2.51 

1157.23 

18.00 

1157.200 

GROSS    MACHINE    UTILIZATION    =       93.61 
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r-l 

/iVPs^AGE    N'H.     np     ITEN.S     PI    THE    QUEUE    WAS 
MAXIMUM  5 


0.339 


QUEUE  CONTENTS 


THE  QUEUE  IS  EMPTY 


CLEUE  PRINTOUT,  QUEUE  NO.   6 

AVERAGE  NO.  OF  ITE^'S  IN  THE  QUEUE  WAS       0.433 
MAXIMUM        5 


QUEUE  CONTENTS 


THE  QUEUE  IS  EMPTY 


CLEUE  PRINTOUT,  QUEUE  NO.   7 


AVERAGf^  NO.  OP  ITEMS  IN  THE  QUEUE  WAS 
MAXIMUM        5 


0.367 


QUEUE  CONTENTS 


THE  QUEUE  IS  EMPTY 


CLEUE  PRINTOUT,  QUEUE  NO.   8 


AVERAGE  NO.  OF  ITEMS  IN  THE  QUEUE  WAS 
NAXIMUM        5 


0.407 


QUEUE  CONTENTS 


THE  QUEUE  IS  EMPTY 


CLEUE  PRINTOUT,  QUEUE  NO.   9 


AVERAGE  NO.  HP  ITlMS  IN  THE  CUEUE  WAS 
NAXIKUV         6 


0.314 


QUEUE  CONTENTS 


THE  QUEUE  IS  EMPTY 


CLEUE  PRINTOUT,  QUEUE  NO.  10 


AVERAGE  NO.  OF  ITFKS  IN  THE  QUEUE  WAS 
^'AXIMU^^        5 


0.320 


QUEUE  CONTENTS 


THE  QUEUE  IS  EMPTY 


CLEUE  PRINTOUT,  GUEUE  NO.  il 


AVERAGE  NO.  OF  ITFN-S  IN  THE  QUEUE  WAS 
^AXIMUM       12  .' 


4.499 


QUEUE  CONTENTS 


1147.47 
1156.99 
1156.99 
1157.49 
1157.49 


1.00  1147.47 
1.00  1156.99 
1.00  1156.99 
l.CO  1157.49 
1.30  1157.49 


1.0€ 

11.00 

0.40 

1297.00 

l.OD 

7.00 

0.82 

1306.00 

l.CO 

7. CO 

1.40 

1306,00 

1.00 

12.00 

0.56 

1307.00 

1.00 

12.00 

1.25 

1307.00 

GLEUE  PRINTOUT,  QUEUE  NO.  12 


AVERAGE  NO.  OF  ITEMS  IN  THE  QUEUE  WAS 
MAXIMUM       13 


4.273 


QUEUE  CONTENTS 


1151.39 
1152.89 
1157.49 


1.00  1151.39 
1.00  1152.89 
l.CC  1157.4  9 


2.00 

12.00 

1.04 

1302.00 

2. CO 

8.00 

0.40 

13U3.00 

2. CO 

12.00 

0.45 

1307.00 

CLEUE  PRINTOUT,  QUEUE  NO.  13 


AVERAGE  NO.  OF  ITEMS  IN  THE  GLEUE  WAS 
NAXIMUM       1/^ 


4.384 
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11.00 

O.^j 

1214.62 

1214.22 

1214.22 

0.0 

2.00 

5.35 

1212.98 

1207.63 

201.27 

7.00 

5.00 

2.20 

1212.98 

1210.78 

433.59 

8.00 

9.00 

2.61 

1226.09 

1223.47 

305.86 

5.00 

12.00 

1.25 

1226.09 

1224.83 

1224.83 

0.0 

10.00  3.51    1240.23 

6.00  2.44    1204.91 

8.00  3.07    1226.09 


1236.72 

412.23 

7.00 

1202.46 

400.82 

8.00 

1223.02 

244.60 

1.00 

QuEub  CONTC-NrS 


115  7.A9 
1157,49 
1157. A9 


l.CC  1157.49 
l.C:  1157.49 
l.Ci  1157.49 


3  .  C  V 

12.0'j 

0.51  1307.C'. 

3.0J 

12.00 

1.33  13o7.0 

3  •  Co 

12.00 

0.4u  1307. OC 

GLEUE  PRINTHUT,  GUeUF  MO.  lA 

AVERAGE  NO.  CF  ITEMS  IN  THE  GUEUE  WAS 
MAXIMUM       13 


4.195 


QUEUE    CONTENTS 
1152.89  1.0       1152.89  4. CO  8.00  1.4f    1303.00 


GLEUE    PRINTOUT,    GUEUE    NC.     15 

AVERAGF    NO.    OF     ITFf'iS     IN    THE    QUEUE    WAS 
NAXIMUM  13 


4.303 


QUEUE 

CONTENTS 

1150.39 

1.0-- 

1150.39 

5  .  00 

7.0Q 

1151.39 

1.00 

1151.39 

5.00 

12.00 

1156.99 

l.OC 

1156.99 

5.00 

7. CO 

1157.49 

i.or 

1157.49 

5.0v 

12.00 

0.43  1300.00 
1.06  1302.00 
0.93  1306.00 
0.40  1307.00 


GLEUE  PRINTOUT,  QUEUE  NO.  16 

AVERAGE  NO.  OF  ITEMS  IN  THE  QUEUE  WAS 
^'AXlMUM       12 


4.450 


QUEUE  CONTENTS 


1156.99 
1157.49 


1.0-  1156.99 
l.CL  1157.49 


6. 00 
6  .  C  0 


7.0G 
12.00 


1306.00  [ 


0.40 

1.01  1307.00 


GLEUE  PRINTOUT,  QUEUE  NO.  17 

AVERAGE  NO.  OF  ITEMS  IN  THE  QUEUE  WAS 
MAXIMUM       13 


4.172 


1151.39 
1156.99 
1157.49 


l.CO  1151,39 
l.CO  1156.99 
l.O:  1157,49 


QUEUE  CONTENTS 

7. GO  12.00 
7. CO  7.00 
7.03    12.00 


1.40  1302.00 

1.40  1306. re 

0.40  13u7.C. 


GLEUE  PRINTOUT,  GUEUE  NO.  18 

AVERAGF  NO,  OF  ITEMS  IN  THE  GUEUE  WAS 


4.253 


126 


3.00 

6.82 

1226.09 

1219.27 

121.92 

6.0  J 

6.0C 

4.89 

1226.09 

1221. 2w 

174.45 

9.00 

1.00 

1.65 

1226.09 

1224.43 

612.21 

1.00 

8.00 


1.4U    1204.91    1203.51    1203.51 


0.0 


6.00 

1.83 

1209.32 

1207.49 

603.74 

8.00 

12.00 

1.06 

1240.23 

1239.16 

1239.16 

0.0 

4.00 

3.13 

1212.98 

1209.85 

302.46 

1.00 

10.00 

2.05 

1226.09 

1224.03 

408.01 

3.00 

7.00     0.40  1212.98  1212.58  1212.58     0.0 
4.00     6.30  1226.09  1219.78   135.53     9.00 


11.00  2.46  1240.23 
3.00  4.53  1212.98 
2.00     7.22  1226.09 


1237.76 

618.88 

5.00 

12C8.45 

241.69 

5.00 

1218.87 

110.80 

3.00 
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**  KFSULTS  OF  THE  STHULATION  RUN  ^-^ 


PRIORITY  RULE  CODE  1 

TOTAL  NO.  OF  JOBS  ARRIVED  1307 
TOTAL  NO.  OF  JOBS  CGMPLETEC  1300 
TOTAL  NO.  OF  JOBS  LATE  0 

TOTAL  UNCOMPLETED  JOBS  IN  SHOP    7. 
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APPENDIX   C 

COMPUTER  PROGRi\M  FT.,OW   CHARTS   AND   LIS  t.INGS 

This    appendl::   includes    the   flow   charts    of   tha   coi?.puter  progriim 
(external  subroutines).      The   coiuputer  progran  listings   with    the 
defiixitinns   of   externcil  subroutines    and  KON-GASP  variables   are  also 
iicludod. 
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flAIN  PROGRAM 


I 


Initialize  Ij'on-G^.ST)  Vcria.bles 


Cell  GASP 


Call  EXIT 


End 


Fig'ore  c-l<  Main  Prograni 
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L 


New  Arrival 

.„ W__,^ 

Gall  ARRVr. 


r 


i 


Return 


End  of  Processing 


Figure  C-2.  Subroutine  EVENTS 


r" 


!    JBFTI.E 


X 


Increase  job  counter 
l)j  one 


T 


Generate  arrival  time  of  the  job 
from  the   siDecified  distribution 


i 


Generate  number  of  operations  and 
the  processing  times  from 
the  specified  distribution 
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initialize  storage  array  STORE 


Generate  the  machine  ordering 
from  the  specified  distribution 


Compute  job  due- date ,  job  slack, 
and  slack  par  remaining  operations 


t 


Store  job   attributes  in 
storage  s^ray  STORE 


return 


F-;r:urc:C-3.  Sub  routine-  .TBFILE 
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ARRYL 


X 


\es 


^ 


Collect  statistics 
for  the  job 


I 


Yes 


w 


Collect  statistics 
fcr  the  machine 


i 


j       Schedule  job  on  the 
I       available  Piachine 


Y.^„„ 


r 


Shop  in  steady 
state  condition? 


> 


i'JO 


Machine  available 
in  specified  machine 
group? 


*v7 


Call  JBFIL 


No 


\V 


Place  job  in 
the  proper  queue 


J 


ope:'ations   of  the  new  job 
in  the  proper  fi  ' -r^- 


Return 


i''i;-;u: 


^+.  Subroutine  ARI^VL 
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U 


'ROC 


TIT 


Identify    . '..'.;   i.iisber 
and  its   'attributes 


I 


Decrease  operations 
for  the  job  by  one 


^es 


1 


I   Collect  statistics 


for  the  job 


t 


Last  operation 
of  the  lob? 


Shop  in  steady 
state  conditions? 


Job  saicple 
completed? 


Print  statistics   / 
collected     / 


N 

\        Return    j 
^^^ ^^ 


No 


i  2    ] 


Figure  C-5.  Subroutine  EPFOC 
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ies 


r 
1.. 


Schedule  Job  on 
the  machine 


©■ 


O 


Find  the  next  operation 

for  the  job 


Machine   available 
in  the  specified 
rnachine  group? 


-£5> 


Jobs  waiting 
in  the  aueue? 


Schedule  job  on 
the  machine 


"1 


1 


Ho 


1 


Place  job  in  Uie 
proter  queue 


No 


Collect  statistics   | 
for  the  i-iac"^_j  r.e    \ 


i 


(    Return    ? 

V        J 


Place  job  operation 
in  the  event  file 


(       n-3tarn        "i 


Figure  C- 5.  Subroutine  EPKOC   (contd.) 
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RULE 

Remove  first  er.bry  from 
the  file  representing  the 
queue  of  the  xaachine 


f    Return    ) 


Figure  C- 6.  Subroutine  RULE 


OUTPUT 


Print  code  for 
the  scheduling  rule 


V 


Print  niUTibers  of  ^obs   arrived, 
jobs   corr.pleted,   jobs   late,   and 
jobs  remaining  in  tVie  shop 


J 


1 


"-s^ 

(        Return         ] 
V_„_ ^ 


Figux'e  C-7.  Subroutine  OITITIIT 
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C  COMPUTER  PROGRAM  FOR  TESTING  SCHEDULING  RULES  IN  AN 

C  INDEPENDENT  NETWORK  OF  WORK  CENTERS  IN  A  JOB  SHOP. 

C  THIS  SIMULATION  PROGRAM  CONSISTS  OF  TWO  GROUPS  OF 

C  ROUTINESf  THE  EXTERNAL  PROGRAM  AND  GASP  ASSEMBLY  PROGRAM, 

C 

C 

C  EXTERNAL  PROGRAM  SUBROUTINES 

C 

C  MAIN  :    THIS  ROUTINE  INITIALIZES  THE  CONTROLLED 

C  VARIABLES  AND  CALLS  SUBROUTINE  GASP. 

C 

c 

C  EVENTS  :  THIS  SUBROUTINE  IS  CALLED  BY  GASP  AND  FORMS  A 

C  LINK  BETWEEN  EXTERNAL  SUBROUTINES  AND  GASP 

C  DECK  .  A  COMPUTED  GO  TO  STATEMENT  CALLS  IN 

C  PROPER  EVENT  SUBROUTINE  DEPENDING  UPON  THE 

C  NATURE  OF  EVENT  . 

C 

C  ARRVL  :   THIS  SUBROUTINE  IS  CALLED  BY  EVENTS  WHENEVER 

C  A  NEW  JOB  ARRIVES  IN  THE  SHOP  .  A  CUSTOMER 

C  WAITS  IN  PROPER  LINE  OR  IS  SCHEDULED  ON  THE 

C  MACHINE  DEPENDING  UPON  WHETHER  THE  MACHINE  IS 

C  BUSY  OR  IDLE. THIS  SUBROUTINE  ALSO  CALLS  JBFILE 

C  WHICH  GENERATES  A  FUTURE  JOB  WHICH  IS  TO  ARRIVE 

C  LATTAR  IN  THE  SHOP. 

C 

C  EPROC  :   THIS  SUBROUTINE  IS  CALLED  BY  SUBROUTINE  EVENTS 

C  WHENEVER  AN  END  OF  SERVICE  EVENT  TAKES  PLACE  . 

C  IT  PROVIDES  INFORMATION  AS  REGARDS  TO  THE 

C  MACHINE  WHICH  IS  TO  PROCESS  NEXT  THE  JOB  AS  PER 

C  THE  TECHNOLOGICAL  ROUTING  SPECIFIED  .  IT  ALSO 

C  CALLS  SUBROUTINE  RULE  WHICH,  DEPENDING  UPON 

C  THE  PRIORITY  RULE,  SELECTS  A  JOB  FROM  A  NUMBER 

C  OF  JOBS  WAITING  IN  A  LINE  OF  THE  MACHINE  UNDER 

C  CONSIDERATION. 

C 

C  RULE  :    THIS  SUBROUTINE  IS  CALLED  BY  EPROC  WHENEVER  A 

C  A  JOB  IS  TO  BE  SELECTED  FROM  A  WAITING  LINE  OF 

C  A  MACHINE  AS  PER  THE  PRIORITY  RULE  SPECIFIED  . 

C 

C  JBFILE  :  THIS  SUBROUTINE  IS  CALLED  WHEN  A  NEW  JOB 

C  ARRIVES  IN  THE  SHOP  .  IT  GENERATES  RANDOMLY 

C  FROM  THE  DESIRED  DISTRIBUTIONS,  VARIOOUS 

C  ATTRIBUTES  OF  THE  JOB  WHICH  ARE  STORED  IN  AN 

C  ARRAY  STORE. 

C 

C 

C  ATTRIBUTES  OF  A  GENERATED  JOB: 

C 

C  ATTRIB(l)     SCHEDULED  TIME  OF  NEXT  EVENT 
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C  ATTRIB(2)     EVENT  CODE 

C  1  =  ARRIVAL 

C  2  =  END  OF  PROCESSING 

C 

C  ATTRIBO)      ARRIVAL  TIME  OF  JOB  IN  SHOP 

C  ATTRIB(^)      MACHINE  GROUP  NUMBER 

C  ATTRIB(5)     NUMBER  OF  OPERATIONS  ON  JOB 

C  ATTRIB(6)     PROCESSING  TIME  OF  OPERATION 

C  ATTRIB{7)     JOB  NUMBER 

C  ATTRIB(8)     OPERATION  NUMBER 

C  ATTRIB(9I     WORK  REMAINING  OF  JOB  ON  ROUTING 

C  ATTRIB(IO)    DUE  DATE  OF  JOB 

C  ATTRIB(ll)    NEXT  MACHINE  fO  OPERATE  ON  JOB 

C 

C 

C   CONTROL  CARDS  VARIABLES  READ  IN  SUBROUTINE  DATAIN: 

C 

C     CONTROL  CARD  NO.  1 

C  NAMF      NAME  OF  THE  PROGRAMMER 

C  NPROJ     PROJECT  NUMBER 

C  MON       MONTH 

C  MDAY      DAY 

C  NYR       YEAR 

C  NRUNS     NUMBER  OF  SIMULATION  RUNS 

C 

C     CONTROL  CARD  NO.  2 

C  NPRAMS    NUMBER  OF  SETS  OF  PARAMETERS  TO  BE  READ  IN 

C  (  NPRAMS  <  20  ) 

C  NHISTO    NUMBER  OF  HISTOGRAMS  TO  BE  OBTAINED  IN 

C  SIMULATION  •  (  NHISTO  <  5  ) 

C  NCOLCT    NUMBER  OF  SETS  OF  STATISTICS  TO  BE  COLLECTED 

C  IN  SUBROUTINE  COLECT  .  (  NCOLCT  <  10  ) 

C  NSTAT     NUMBER  OF  SETS  OF  STATISTICS  TO  BE  COLLECTED 

C  IN  SUBROUTINE  TMSTAT  .  (  NSTAT  <  10  ) 

C  ID        NUMBER  OF  COLUMNS  IN  ARRAY  NSET 

C  IM        NUMBER  OF  ATTRIBUTES  ROW  IN  NSET 

C  NOO       NUMBER  OF  FILES  TO  BE  MAINTAINED 

C  MXC       LARGEST  NUMBER  OF  CELLS  USED  IN  HISTOGRAM 

C  (  MXC  <  22  I 

C  SCALE     A  PARAMETER  USED  TO  MULTIPLY  FLOATING  POINT 

C  ATTRIBUTES  PRIOR  TO  CHANGING  THEM  TO  FIXED 

C  POINT  FOR  STORAGE  IN  NSET 

C 

C     CONTROL  CARD  NO.  3 

C  MSTOP     CODE  FOR  METHOD  OF  ENDING  SIMULATION 

C  MSTOP  =  OtSIMULATION  TERMINATES  WHEN  END 

C  OF  SIMULATION  EVENT  OCCURS 

C  MSTOP  >  0, SIMULATION  TERMINATES  WHEN 

C  PREDETERMINED  TIME  IS  REACHED 

C  JCLEAR    JCLEAR  -    0  t  SIMULATION  IS  STARTED  FROM 
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C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


BEGINN 

CHANGI 

JCLEAR 

PRIOR 

NORPT 

NORPT 

ARE  BY 

NORPT 

NEP 

NEP  = 

DATAIN 

NEP  = 

ARRAYS 

TSTART 

STARTI 

TSTOP 

TIME  T 

NSEED 

NUMBER 

AND  WITHOUT 


ARE  USED 
SUBROUTINE 

STATISTICAL 


ING    (  TNOW  =  TSTART  ) 
NG  COMDITIONS  OF  RUN 

=  1  .  SYSTEM  IS  INITIALIZED  AGAIN 
TO  REPEATING  A  SIMULATION  RUN 
=  1  •  SUBROUTINES  SUMARY  &  OUTPUT 
PASSED 

=  0  f  SUMARY  &  OUTPUT 
1  ,  NEW  PARAMETERS  IN 

ARE  READ  IN 
1  ,  ONLY  NSET,  TNOW  G 

ARE  INITIALIZED 
NG  TIME  OF  SIMULATION 
0  END  SIMULATION 

OF  SEED  VALUES  TO  BE  READ  IN 


CONTROL  CARD  NO,  4 

IX(I)     VALUE  OF  INITIAL  SEED  NUMBER  USING  SEED 
NO.  I  FOR  GENERATING  JOB  ATTRIBUTES 

CONTROL  CARD  NO.  5 

INN(I)        IF  INNU)  =  1,  THE  ENTRIES  IN  FILE  I  ARE 
ORDERED  BY  KRANKd)  FROM  LOWEST  VALUE  TO 
HIGHEST  VALUE(FIFO),  IF  INN(n  =  2  8  THE 
ENTRIES  ARE  ORDERED  BY  KRANKd)  FROM 
HIGHEST  VALUE  TO  LOWEST  VALUE(LIF0)9 


CONTROL  CARD  NO.  6 

KRANKd)  JOB  ATTRIBUTE  ROW  ON  WHICH  FILE 


I  IS  RANKED 


CONTROL 
JSMP 
MGRP 
NRUL 


KODE 


NJOB 


CARD  NO.  7 

JOBS  PER  SIMULATION  RUN 

NO.  OF  MACHINE  GROUPS  IN  SHOP 

PRIORITY  RULE  CODE  NUMBER 

1  SHORTEST  IMMINENT  OPERATION 

2  FIRST  COME  FIRST  SERVED 

3  EARLIEST  JOB  DUE  DATE 

4  JOB  SLACK  TIME 

5  SLACK  PER  REMAINING  OPERATIONS 

CODE  FOR  METHOD  OF  COMPUTING  DUE  DATES 
2   MULTIPLIER  EQUAL  TO  NUMBER  OF 
OPERATIONS  PER  JOB 

JOB  ARRIVAL  COUNTER 


CONTROL  CARD  NO.  8 

MACH(I)   NUMBER  OF  MACHINES  IN  MACHINE  GROUP  I 

CONTROL  CARD  NO.  9 

NCELLStI)  NUMBER  OF  CELLS  IN  HISTOGRAM  I 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 


CONTROL  CARDS  NO.  10,11,12,13 

EACH  OF  THE  FOLLOWING  CARD  CONSISTS  OF  FOUR 
PARAMETERS,  J=l,2,3,4 

-  ARRIVAL  DISTRIBUTION  PARAMETERS  (1=1) 

-  INTERARRIVAL  DISTRIBUTION  PARAMETERS  (1=2) 

-  NUMBER  OF  OPERATIONS  DISTRIBUTION  PARAMETERS  (1=3) 

-  MACHINE  ROUTING  DISTRIBUTION  PARAMETERS  (1=4) 

-  PROCESSING  TIME  DISTRIBUTION  PARAMETERS  (1=5) 


CONTROL  CARD  NO.  11 


KDSTAR 
KDSTIA 
KDSTOP 
KDSTRT 
KDSTPR 


CODE 
CODE 
CODE 
CODE 
CODE 


FOR 
FOR 
FOR 
FOR 
FOR 

1 
2 
3 

5 


ARRIVAL  DISTRIBUTION 
INTERARRIVAL  DISTRIBUTION 
DISTR.  FOR  NUMBER  OF  OPERATIOS 
MACHINE  ROUTING 
PROCESSING  TIME  DISTRIBUTION 

ERLANG 

POISSON 

GEOMETRIC 

UNIFORM 

NORMAL 


CONTROL  CARDS  NO.  14 

THE  POSITION  OF  THE  CONTROL  CARD(  WITHIN  CONTROL 
CARDS  NO.  10,11,12,13  )  TO  BE  READ  FOR  SPECIFYING; 
ICDARR    ARRIVAL  DISTRIBUTION 
ICDIAR    INTERARRIVAL  DISTRIBUTION 
ICDOPR    NUMBER  OF  OPERATION  DISTRIBUTION 
ICDRTG    MACHINE  ROUTING  DISTRIBUTION 
ICDPRC    PROCESSING  TIME  DISTRIBUTION 


EXTERNAL  PROGRAM  VARIABLES 


LTNES 

TBD 

JA 

JD 

MG 

lOPN 

JLAT& 

NOPRM(J) 

MACH(I) 

MAVL(I) 

MCG( J) 

XIAT 

TLA 


JOB  LATENESS 
TIME  BETWEEN  DEPARTURES 
JOB  ARRIVAL  COUNTER 
JOB  DEPARTURE  COUNTER 
NUMBER  OF  MACHINE  GROUPS 

OPERATIONS 

JOBS  LATE 

OPERATIONS  REMAINING  FOR  JOB 

MACHINES  IN  MACHINE  GROUP  I 


IN  THE  SHOP 


NUMBER 

NUMBE- 

NUMBER 

NUMBER 

NUMBER 

GROUP  I 

MACHINE  NUMBER  FOR 

INTER-ARRIVAL  TIME 

TIME  OF  LAST  ARRIVAL 


OF 
OF 
OF 
OF 
OF 


MACHINES  AVAILABLE  IN  MACHINE 


JOB  J 


OF  A  JOB 
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C  DDT  OUF-DATE  OF  A  JOB 

CC  XISHP  NUMBER  OF  JOBS  IN  THE  SHOP 

C  XMBUZ  NUMBER  OF  MACHINES  BUSY 

C  TLD  TIME  OF  LAST  DEPARTURE  OF  A  JOB 

C  TBD  TIME  BETWEEN  DEPARTURES  OF  JOBS 

C  SPTM(J)  SUM  OF  PROCESSING  TIMES  OF  JOB  J 

C  SHOP(J)  SHOP  FLOW  TIME  OF  JOB  J 

C  PTM(J)  PROCESSING  TIME  OF  JOB  J 

C  MGRP  NUMBER  OF  MACHINE  GROUPS  IN  THE  SHOP 

C  NPRL  PRIORITY  RULE  NUMBER 

C 

c 
c 

C  MAIN  PROGRAM 


C 


COMMON  ID,IM,INIT,JEVENT, JMONIT,MFA,MSTOP,MX,MXC, 
1N0RPT,N0T,NPRAMS,NRUN,NRUNS,NSTAT,NSEED,MXX,NEP, 
2NC0LCT,NHIST0,N0Q 

COMMON  INN(19), JCELLS(5,22),KRANK(I9),MAXN0(19) , 
1MFE(19),MLC(19),MLE(19),NCELLS(5),N0( 19),IX(8) 

COMMON  OUT,SCALE,TNOW,TSTART,TSTOP 

COMMON  ATTRIB(13),ENQ(19),PARAMS(20,4),QTIME(19>, 
1SSUMA(3C,5),SUMA(50,5) 

COMMON  MGQ,LTNES,NJOB,JA, JD, JOB,KOOE,KTR,MG, 
IJWRM, JSMP, ISUMR,IOPN,LONE,  JLATE 

COMMON  N0PRM(150C  )  ,MACHa9  )  ,MAVL  ( 19  ) ,  MCG(  20  ) 

COMMON  WRM, XI  AT, TLA, DDT, XISHP, XMBUZ, TLD, TBU 

COMMON  SPTM{1500) , STORE ( 15 , 15 ), TSHOP { 1500 ) , 
1W0RK0{15),PTM(25) 

COMMON  MGRP, NPRL, NRUL, ICDARR, ICDIAR, ICOOPR, 
1ICDRTG,ICDPRC,KDSTAR,KDSTIA,KDST0P,KDSTRT,K0STPR 

COMMON  XLTNS 

DIMENSION  NSET(15,5000) 
JLATE  =  0 
JD  =  0 
TLA  =  0. 
TLD  =  0. 
ATTRIB(l)  =  0. 
BT  =  0. 
XISHP  =  1. 
DO  20  1=1,1500 
NOPRM(I)  =  0 
TSHOP( I)  =  0, 
20  CONTINUE 

CALL  GASP(NSET) 

CALL  EXIT 

END 


142 


SUBROUTINE  JBFIL 

c 

c 

ERLANG       1 

c 

POISSON      2 

c 

GEHMETRIC    3 

c 

UNIFORM      4 

c 

NORMAL       5 

c 

c 

COMMON  ID, IM,INIT, JEVENT,JMONIT,MFA,MSTOP,MX»MXC, 
lNORPT,NOT,NPRAMS,NRUN,NRUNS,NSrAT,NSEED,MXX,NEP, 
2NC0LCT,NHIST0,N0Q 

COMMON  INN(19),JCELLS(5,22) ,KRANK(I9) ,MAXNQ(19), 
1MFE{I9),MLC(19),MLE(I9),NCELLS(5),NQ(19),IX{8) 

COMMON  OUT,  SCALE,  TNOVy,TSTART,TSTOP 

COMMON  ATTR IB ( 13 ),ENO( 19), PAR AMS{ 20,4 ),QTIMEi 19), 
1SSUMA(30,5) , SUM A (50,5) 

COMMON  MGQ,LTNES,NJOB,JA,JD,JOB,KOOE,KTR,MG, 
IJWRM, JSMP, ISUMR,IOPN,LONE, JLATE 

COMMON  NOPRM(150D),MACH(i9),MAVL(19),MCG(20) 

COMMON  WRM,XIAT,TLA,DDT,XISHP,XMBUZ,TLD,TBO 

COMMON  SPTM( 1500), STORE ( 15 , 15 ) , TSHOP ( 1500 ) • 
1W0RK0(15) ,PTM(25) 

COMMON  MGRP,NPRL,NRUL,ICDARR, ICOI AR, ICOOPR, 
1ICDRTG,ICDPRC,KDSTAR,KDSTIA,KDST0P,K0STRT,KDSTPR 

COMMON  XLTNS 

DIMENSION  NSET(15,1) 
ATTRIB(l)  =  0. 
NJOB  =  NJOB+1 
JA  =  NJOB 

IF(KOSrAR.EO.O)  GO  TO  5 
GO  TO  (1,2,3),KDSTAR 

1  PRINT  4000 

4000  F0RMAT(1H1,«  WRONG  DIST.  SELECTED  •/) 
CALL  EXIT 

2  ARVTM  =  NP0ISN(1,IC0ARR) 
GO  TO  14 

3  ARVTM  =  GEMTR(1,ICDARR) 
14  ATTRIB(l)  =  ARVTM 

GO  TO  16 

5  IF(KDSTIA.EQ.O)  GO  TO  6 
GO  TO  7 

6  PRINT  406 

406  FORMATdHl,'  NO  DIST.  SELECTED  FOR  ARRIVALS  •/) 
CALL  EXIT 

7  GO  TO  (8,1,1,1,9),KDSTIA 

8  BARVTM  =  ERLANGd,  ICDIAR) 
GO  TO  10 

9  BARVTM  =  RNORMLCl, ICDIAR) 


143 


10    GO    TO    12 

12    ATTRIB(l)     =    BT-i-BARVTM 
BT    =    ATTRIB(l) 

16  ATTRIB(2)     =    1. 
ATTRIBO)     =    ATTRIB(l) 
ATTRIB(7)     =    JA 

C 

C     GENERATING  NUMBER  OF  OPERATIOS  PER  JOB 

C 

GO  TO  {1,1,18,17,1),KDST0P 

17  A  =  PARAMS( ICDOPR.l) 
B  =  PARAMS(ICD0PR,2) 
lOPN  =  A+(B-A)*DRAND(1) 
GO  TO  19 

18  lOPN  =  GEMTR( 1, ICDOPR) 

19  ATTRIB(5)  =  lOPN 
C 

C     GENERATING  PROCESSING  TIME  OF  EACH  OPERATION 
C 

SPTM(JA)  =  0, 

GO  TO  (21,1,1,1,22)»KDSTPR 

21  DO  20  I=1,I0PN 

PTMII)  =  ERLANGCl, ICDPRC) 
SPTM(JA)  =  SPTM(JA)+PTM( I) 
DO  2C  J=l,ll 
STORE! I, J)  =  0. 

20  CONTINUE 
GO  TO  24 

22  DO  23  1=  1,I0PN 

PTM(I)  =  RN0RML(1,ICDPRC) 
SPTM(JA)  =  SPTM{JA)+PTM( I ) 
DO  23  J=l,ll 
STORE! I, J)  =  0. 

23  CONTINUE 
C 

C      GENERATING  MACHINE  ROUTING 
C 

24  WRM  =  SPTM(JA) 
IP  =  10 

DO  60  1=1, lOPN 
35  MCG(I)  =  1+(MGRP)*DRAND( 1) 
IF(MCG(I) .EO.IP)  GO  TO  35 
ATTRIB(A)  =  MCG(I) 
IP  =  MCG( I ) 
ATTRIB(8)  =  I 
ATTRIB(6)  =  PTM(I) 
IFd.EO.l)  GO  rO  40 
WRM  =  WRM-PTM( I-l) 
ATTRIB(9)  =  WRM 
GO  TO  (36,37) ,KODE 
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36  ATTRIB(IO)  =  ATTRIB ( 1 )*9.*SPrM( JA) 
GO  TO  50 

37  ATTRIB(IO)  =  ATTR  I  B  (  1  )-HOPN«SPTM  (  JA) 
GO  TO  50 

40  ATTRIB{9)  =  WRM 

GO  TO  (38,39) ,KODE 

38  ATTRIB(IO)  =  AT  TR  I  B  ( 1 )  •i-9.*SPTM(  J  A) 
GO  TO  5C 

39  ATTRIB(IO)  = 
50  ATTRIB(ll)  = 

ATTRIB(12)  = 

DO  60  J=1,12 
60  STORE(  I, J)  = 

STORFt I0PN,13)  = 

LL  =  IOPN-1 

DO  70  1=1, LL 
70  STORE(  1,13)  =  MCG( I  +  l) 

RETURN 

END 


ATTRIB( I)^IOPN*SPTM( JA) 

ATTRIB(10)-ATTR£8(9) 

ATTRIB(11)/(I0PN-I+1) 

ATTRIBtJ) 

0. 
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SUBROUTINE  EVENTS ( I »NSET ) 

COMMON  ID, IM, INIT,JEVENT, JM0NIT,MFA,MST0P,MX,MXC, 
lNORPT,NOT,NPRAMS»NRUN,NRUNS,NSTAT,NSEEDtMXX,NEP, 
2NC0LCT,NHIST0,N00 

COMMON  INN(19),JCfcLLS(5,22),KRANK(19),MAXNQ(I9), 
1MFE(19),MLC(19),MLE(I9),NCELLS(5»,NQ(19),IXC8) 

COMMON  OUT, SCALE, TNOW,TST ART, TSTOP 

COMMON  AT TR IB (13) ,ENQ (19 ) , PAKAMS ( 20,4 ) , OT I  ME ( 19 ) , 
1SSUMA(30,5),SUMA(50,5) 

COMMON  MGQ,LTNES,NJOB,JA,JD,JOB,KODE,KTR,MG, 
IJWRM, JSMP,ISUMR, I OPN, LONE , JLATE 

COMMON  NOPRM( 1500),MACH( 19) ,MAVL(19) ,MCG(20) 

COMMON  WRM,XIAT,TLA,DDT,XISHP,XMBUZ,TLO,TBD 

COMMON  SPTM ( 1 500 ), STORE ( 15,1 5 ),TSHOP (1500), 
1W0RKQ(15),PTM(25) 

COMMON  MGRP,NPRL,NRUL,ICDARR,ICOIAR,ICDOPR, 
1ICDRTG,ICDPRC,KDSTAR,KDSTIA,KDST0P,K0STRT,KDSTPR 

COMMON  XLTNS 

DIMENSION  NSET(15,1) 
GO  TO  (1,2), I 
CALL  ARRVL(NSET) 
RETURN 
:  CALL  EPROC(NSET) 
RETURN 
END 
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SUBROUTINE  ARRVL(NSET) 
C 

COMMON  ID, IM, INIT,JEVENT,JMONIT,MFA,MSTOP,MX,MXC, 
1N0RPT,N0T,NPRAMS,NRUN,NRUNS»NSTAT,NSEED,MXX,NEP, 
2NC0LCT,NHIST0,N00 

COMMON  INN(19),JCELLS(5,22) ,KRANK( 19) tMAXNOi 19 ) , 
1MFE(19),MLC(19) ,MLE ( 19 ) , NCELLS ( 5 ) , N0( 19 ) , I  X ( 8 ) 

COMMON  OUT,SCALE,TNOW,TSTART,TSTOP 

COMMON  ATTRIB(13),ENQ(19),PARAMS«20,4),QTIME(19), 
1SSUMA(30»5),SUMA(50,5) 

COMMON  MG0,LTNES,NJ08, JA, JD,JOB,KODE,KTR,MG, 
IJWRM, JSMP, ISUMR,IOPN,LONE, JLATE 

COMMON  NOPRM ( 1500 ) , MACH( 19 ) , MAVL ( 19 ) , MCG( 2C ) 

COMMON  WRM,XIAT,TLA,ODT,XISHP,XMBUZ,TLD,TBD 

COMMON  SPTMI 1 500 ), STORE (15,1 5 ),TSHOP (1500), 
1W0RKQ(15) ,PTM(25) 

COMMON  MGRP,NPRL,NRUL, ICDARR, ICDIAR, ICDOPR, 
1ICDRTG,ICDPRC,KDSTAR,KDSTIA,KDST0P,KDSTKT,K0STPR 

COMMON  XLTNS 
C 

DIMENSION  NSET( 15,1) 

CALL  TMSTAT(XISHP,TN0W,10,NSET) 

XISHP  =  XISHP+1. 

XIAT  =  TNOW-TLA 

CALL  HISTOG(XIAT,0.,0.4,1) 

CALL  C0LECT(XIAT,42,NSET) 

TLA  =  TNOW 

LSF  =  1 

DO  20  1=101,1501,100 

LSF  =  LSF^-1 

IF( JD.LT. I )  GO  TO  21 

20  CONTINUE 

21  IF( JD.GT.800)  GO  TO  23 

22  JSF  =  LSF+Al 

CALL  C0LECT(XIAT,JSF,NSET) 

23  CALL  HISTOG(TNOW,0., 0.4,2) 
lOPN  =  ATTRIB(5) 
NOPRM(JA)  =  lOPN 

MG  =  ATTRIB(4) 

IF  (MAVL(MG) )12,30,25 
12  CALL  ERROR(NSET) 
C 

C         SCHEDULE  NEXT  JOB  ON  IDLE  MACHINE 
C 

25  XMBUZ  =  MACH(MG)-MAVL(MG) 

CALL  TMSTAT(XMBUZ,TNOW,MG,NSET) 

MAVL(MG)  =  MAVL(MG)-1 

ATTRIB(l)  =  TN0W+ATTRIB(6) 

ATTRIB(2)  =  2, 

CALL  FILEM(1,NSET) 
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GO  TO  65 
C 

C  IF  NO  MACHINE  IS  AVAILABLE  THEN  JOB  MUST  WAIT 

C 

30  ATTRIB(I)  =  TNOW 

ATTRIB(2)  =  1. 

MGQ  =  MG+1 

CALL  FILEM(MGQ,NSET) 
C 

C         RELEASE  NEXT  JOB  IN  THE  SHOP 
C 

65  CALL  JBFILE(NSET) 

DO  70  K=l,13 
70  ATTRIB(K)  =  ST0RE(1,K) 

ATTRIBO)  =  ATTRIB(I) 

lOPN  =  ATTRIB(5) 

CALL  FILEMd.NSET) 

DO  80  I=2t lOPN 

DO  90  J=l,13 
90  ATTRIBt J)  =  STORE( I, J) 

ATTRIB(3)  =  ATTRIB(l) 

MG  =  ATTRIB(4) 

LONE  =  MG-HO 

CALL  FILEM{LONE,NSETI 
80  CONTINUE 

RETURN 

END 
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SUBROUTINE  RULECNSET) 

COMMON  ID,IM,INIT,JEVENT,JMONIT,MFA,MSTOP,MX,MXC, 
lNORPT,NOT,NPRAMS,NRUN,NRUNStNSTAT,NSEED,MXX,NEP, 
2NC0LCT,NHIST0,N0Q 

COMMON  INN(19), JCELLS(5,22) ,KRANK( 19) , MAXNQ ( 19 ) , 
1MFE(19),MLC(19),MLE(19),NCELLS(5),N0(19),IX(8) 

COMMON  OUT,SCALE,TNOW,TSTART,TSTOP 

COMMON  ATTRIB{13),ENQ(19),PARAMS(20,4),QTIME( 19), 
1SSUMA(30,5),SUMA(50,5) 

COMMON  MGO,LTNES,NJOB, JA, JD,JOB,KODE,KTR,MG, 
IJWRM, JSMP, ISUMR,IOPN,LONE,  JLATE 

COMMON  NOPRM(1500),MACH(19),MAVL(19),MCG(20) 

COMMON  WRM,XIAT,TLA,ODT,XISHP,XMBUZ,TLD,TBD 

COMMON  SPTM(1500),ST0REa5,15),TSH0P(1500), 
1W0RKQ( 15) ,PTM(25) 

COMMON  MGRP,NPRL,NRUL, ICDARR, ICDIAR, ICDOPR, 
1IC0RTG,IC0PRC,K0STAR,KDSTIA,KDST0P,KDSTRT,KDSTPR 

COMMON  XLTNS 

DIMENSION  NSET(15,1) 

MFFQ  =  MFF(MGQ) 

CALL  REMOVE(MFEQ,MGQ,NSET) 

RETURN 

END 
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SUBROUTINE  EPROC(NSET) 
C 

COMMON  ID, IM, INIT,JEVENTf JMONIT,MFA,MSTOP,MX,MXC, 
1N0RPT,N0T,NPRAMS,NRUN,NRUNS»NSTAT,NSEEU,MXX,NEP, 
2NCOLCT,NHISTO,NO0 

COMMON  INN(19), JCELLS(5,22),KRANK(19),MAXNQ(19), 
1MFE(19),MLC(19),MLE(19),NCELLS(5),NQ( 19),IX(8) 
COMMON  OUT, SCALE, TNOWt TST ART, TSTOP 

COMMON  ATTRIBC 13) ,ENQ(19) , PARAMS( 20,4 ) , QTI ME ( 19 ) , 
1SSUMA(30,5) ,SUMA(50,5) 

COMMON  MGQ,LTNES,NJOB, JA, JD,JOB,KODE,KTR,MG, 
IJ^iRM,  JSMP,  ISUMR,IOPN,LONE,  JLATE 
COMMON  NOPRM(1500) ,MACH( 19 ) , MAVL ( 19 ) , MCG( 20 ) 
COMMON  WRM,XIAT,TLA,DDT,XISHP,XMBUZ,TLD,TBD 
COMMON  SPTM(1500),STORE(15,15),TSHOP(1500), 
1W0RKQ(15) ,PTM(25) 

COMMON  MGRP,NPRL,NRUL, ICDARR, ICDIAR, ICDOPR, 
1ICDRTG,IC0PRC,KDSTAR,KDSTIA,K0ST0P,K0STRT,KDSTPR 
COMMON  XLTNS 
C 

DIMENSION  NSET(15,1) 
7000  FORMATdH  ,415) 

LONE  =  ATTRIB{13)4-10. 
CDT  =  ATTRIB(IO) 
MG  =  ATTRIB(4) 
MGO  =  MG-H 
JOB  -    ATTRIB(7) 
NOPRM(JOB)  =  N0PRM(J0B)-1 
C 

C  IF  ALL  OPERATIONS  HAVE  BEEN  COMPLETED  ON  A  JOB 

C 

IF  (NOPRM(JOB))  10,10,50 
10  IF  (ATTRIB(3)-TSTART)20,20,5 
5  IF  (ATTRIB(3)-TST0P)2,2,20 
C 

C         COLLECT  STATISTICS  ON  COMPLETED  JOB 
C 

2  M  =  ATTRIB(5) 

XLTNS  =  ATTRIB(1)-DDT 
IF(XLTNS.GT.O)  GO  TO  12 
GO  TO  14 
12  JLATE  -    JLATE+1 

14  CALL  C0LECT(XLTNS,1,NSET) 
CALL  HISTOG{XLTNS,-200.,20.,4| 

15  TSHOP(JOB)  =  ATTRIB(1)-ATTRIB(3) 
CALL  COLECT(TSHOP{ J0B),16,NSET) 
CALL  HISTOG(TSHOP{ J0B),5., 10.,5) 
CALL  TMSTAT(XISHP,TN0W,10,NSET) 
XISHP  =  XISHP-1. 

JD  =  JD+1 
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TBD    =    TNOW-TLO 

CALL    C0LECT(TBD,32,NSET) 

TLD    =    TNOW 

IF( JD.GT.800)    GO    TO    6 
8    XIAVG    =    SUMA(42,1)/SUMA(42,3) 

XBDVG  =  SUMA(32f 1)/SUMA{32,3) 

XTSVG  =  SUMAC 16,1)/SUMA{16,3) 

XLTVG  =  SUMA(1,1)/SUMA{1,3) 
1000  FORMATdH  ,  2  16  ,  F8.  1  ,4F10.2  ) 

PRINT  1000, J A, JO, TNOW, XIAVG, XBDVG, XTSVG, XLTVG 
6  CALL  HIST0G{TBD,0*,0.4,3) 

NSF  =  1 

DO  22  1=101,1501,100 

NSF  =  NSF+l 

IFtJD.LT. I )  GO  TO  23 

22  CONTINUE 

23  MSF  =  NSF+15 

CALL  COLECT(XLTNS,NSF,NSET) 
CALL  COLECT(TSHOP( JOB),MSF,NSET) 
IF( JD,GT,800)  GO  TO  17 
21  KSF  =  NSF^31 

CALL  COLECT(TBD,KSF,NSET) 
GO  TO  17 

16  JD  =  JD+1 
XISHP  =  XISHP-1. 
TBD  =  TNOW-TLD 
TLD  =  TNOW 

17  IF{ JA.EQ.1499)  GO  TO  3 
IF  (JD-JSMP)  20,3,3 

C 

C         JOB  SAMPLE  HAS  BEEN  COMPLETED  AT  THIS  POINT 

C 

3  NEP  =  1 
MSTOP  =-1 
NORPT  =0 
RETURN 
C 

C         UPDATING  STATUS  OF  MACHINE  GROUP 
C 

20  IF  (NOCMGQ))  30,30,40 

30  XMBUZ  =  MACH(MG)-MAVL(MG) 

CALL  TMSTAT(XMBUZ,TNOW,MG,NSET) 

18  MAVL(MG)  =  MAVL(MG)+1 
RETURN 

C 

C         SELECT  NEXT  OPERATION  OF  JOB 

C 

50  XJOB  =  JOB 

CALL  FIND(XJ0B,7,L0NE,5,KC0L,NSET) 
CALL  REHOVE(KCOL,LONE,NSET) 
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MCH  =  ATTRIB(4) 

MCF  =  MCH-H 

IF(MAVL(MCH) )  45,45,60 
45  ATTRIB(l)  =  TNOW 

ATTRIB(2)  =  1. 

CALL  FILEM(MCF,NSET) 
70  IF(NO(MGO))  55,55,40 
55  XMBUZ  =  MACH(MG)-MAVL(MG) 

CALL  TMSTAT(XMBUZ,TNOW,MG,NSET) 
19  MAVL(MG)  =  MAVL{MG)+1 

RETURN 
60  MAVLIMCH)  =  MAVL(MCH)-1 

ATTRIB(l)  =  TN0W+ATTRIB(6) 

ATTRIB(2)  =  2. 

CALL  FILEM(1,NSET) 

GO  TO  7C 
40  CALL  RULE(NSET) 

ATTRIB(l)  =  TN0W*ATTRIB(6) 

ATTRIB{2)  =  2. 

WORKOiMG)  =  WORK0(MG)-ATTRIB(6) 

CALL  FILEM(l,NSeT) 

RETURN 

END 
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SUBROUTINE  OUTPUT(NSET) 


COMMON  ID 
INORPTtNOT 
2NC0LCT,NH 

COMMON  IN 
1MFE(19),M 

COMMON  OU 

COMMON  AT 
1SSUMA( 30» 

COMMON  MG 
IJWRM, JSMP 

COMMON  NO 

COMMON  WR 

COMMON  SP 
1W0RKQ( 15) 

COMMON  MG 
IICDRTG,IC 

COMMON    XL 


, IM, INIT 
,NPKAMS, 
ISTO.NOQ 
N(19), JC 
LC(19),M 
T, SCALE, 
TRIB(13) 
5) ,SUMA( 
Q.LTNES, 
,ISUMR,I 
PRM( I50C 
M,XIAT,T 
TM( 1500) 
,PTM(25) 
RP,NPRL, 
OPRCKDS 
TNS 


,JEVENT,JMONIT,MFA,MSTOP,MX,MXC, 
NRUN,NRUNS,NSTAT,NSEED,MXX,NEP, 

ELLS (5, 22) ,KRANK(19) ,MAXNQ(19) , 

LE(19),NCELLS(5),NQ(19),IX(8) 

TNOW,TSTART,TSTOP 

,EN0(19),PARAMS(20,4),QTIME(19), 

50,5) 

NJ08,JA, JD,JOB,KODE,KTR,MG, 

OPN,LONE, JLATE 

),MACH(19) ,MAVL(19),MCG(20) 

LA,DDT,XISHP,XMBUZ,TLD,TBO 

,STORE(15,15),TSHOP(150G), 

NRUL,ICDARR, ICDI AR, ICDOPR, 
TAR,KDSTIA,KDSTnP,KDSTRT,KDSTPR 


100 
200 
300 
400 
500 
600 
700 


DIMENSION 

FORMATCIH 

FORMAT(  IH 

FORMATdH 

FORMAT (IH 

FORMATdH 

FORMATdH 

FORMATdH 

PRINT  ICO 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

RETURN 

END 
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**  RESULTS  OF  SIMULATION  RUN  ** 
PRIORITY  RULE  CODE  '.yXjIlO/) 
TOTAL  NO.  OF  JOBS  ARRIVED* , 110/ ) 

NO.  OF  JOBS 

NO.  OF  JOBS 


TOTAL 
TOTAL 


COMPLETED  SIT/) 
LATE  •,3X,I10/) 


200 
300 
400 
500 
700 
600 


1,20X,»  **  RESULTS  OF  SIMULATION  RUN  **  •//) 

f20X, 

»20X, 

»20X, 

f20X, 
1) 

,20X,«  TOTAL  UNCOMPLETED  JOBS  IN  SHnP»,F5.0/) 

,NRUL 
,JA 
,  JD 
, JLATE 

,XISHP 
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Tf«e  purpose  of  the  research  ic  to  (1)  present  a  consistent  and 
tnore  co.r.prcaciisive  revitv/  of  the  most  iripcrtant  joh  shop  siiv^ulatlon 
C;xpcrirT:cn'.  s  conducted  since  1960;  (2)  develop  a  method  of  assigning 
job  duG-datcs,  c.nd  testing  it;  (3)  evaluate  the  efficiency  of  GASP-II 
in  bsTidling  .slinulaticn  experiments;  and  (4)  study  the  effect  ("jf  dif~ 
fer'cnt  pfoccr-sing  time  distribu  .ion  on  the  perfonnance  of  the  rules. 

Consid-viable  experimental  investigation  is  conducted.   The 
most  significant  results  are:  (1)  the  shortest  imniiuent  operation 
rule  ii.-.  snperior  to  others  in  reducing  job  lateness  and  shop  flow 
tiTne;  (2)  the  procedure  in  which  the  due-date  allowance  is  proportional 
to  the  nivjhcr  of  operations  and  v;ork  content  of  the  jobs  has  p;cvcd 
to  h2  beneficial  in  the  case  of  the  non-due-date  rules  such  as  first-- 
coiic  f  i  Tst-sevvc  "  and  shortest  imminent  operation  tin;e  rules;  (3)  the 
perforr.ianct  of  li.3  shortest  imjTiinent  operation  tim.e  and  slack  p^er 
remaining  number  jf  operations  rules  are  affected  badly  with  the  use 
of  the  processing  :ime  dis  rributions  of  Erlang  v;ith  parar.ietcr  K  equallir.g 
^  and  8,   However,  the  performances  under  Erlang  8  are  better  than 
those  obtained  ur.der  Erlang  4;  and  (A)  the  GASP-II  vjorks  equally  good 
for  large  size  of  shop  problems.   The  simulation  com.puter  execution 
time  of  runs;  ranging  from  500  to  1300  departed  jobs,  varies  from  18 
to  23  m.inutes  depending  on  the  rules  tested. 


